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Interactive computer systems, forecasting, crisis, crisis and policy 
changes, China 


-This report sunrvirizrs the development and testing of the CEVMKC interactive 
forecasting system for monitoring political crises and policy changes in 
China. The first chapter presents an overview of the project, discusses 
its objective and its major accosipl ishraents and delineates the structure 
of the report. The second chapter discusses the philosophy and the design 
of the CI'WMFC System, and describes the forecasting model and the modeling 
strategy the system employs. The third chapter delineates the structure 






















esults of sensitivity analysis and forecast reliability tests per- 
ormed on the CEWMFC system, mid defines the hypotheses and the 
ndicators employed in data analyses. Appendix A then presents 
technical description of the uvithenatleal foundation of the CEVMPC 


ivstcm. Appendix B describes the coding procedure for constructing 
:he Chinese crisis and policy indicators. Appendix C describes the 
itatistical models used in the data analyses. 
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Introduction 


This report summarizes the results derived from a project to develop 
an interactive computer forecasting system for crisis early warning and 
monitoring. The project was conducted at Michigan State University under 
the sponsorship of the Defense Advanced Research Projects Agency's Cybernetics 
and Technology Office as a part of its Crisis Management Program. 

1.0 Project Overview 

The project was conceived in response to Crisis Management Program’s 
need for auxilllary computer software to augment the forecasting capability 
of its existing Early Warning and Monitoring System (EVIAMS). This 
auxilllary system, designated the Crisis Early Warning and Monitoring 
Forecaster (CEWMFC), generates through leading Indicators multivariate 
forecasts of crisis frequencies and probabilities as distinct from the 
subjective forecasts of the existing system. 

The EWAMS generates, from an inventory of crisis Indicators, estimates 
of the probability that the observed value at a given time point reflects an 
abnormal deviation from "usualness”, which is taken to indicate a crisis 
occurence. Since these probabilities are estimated for the existing 
observations, the analyst must intuit the likelihood of future crisis events. 

Supplementing the foregoing subjective forecasting capability, the CEWMFC 
system offers a frequency forecasting capability through which one could 
numerically extrapolate future crisis probabilities from past trends and 
through coterminous leading indicators. The use of mathematical routines to 
project future values from past trends minimizes the biases attendant to 
subjective evaluation and estimation of future likelihood based on current 


values. 




The CEWMFC system Is also fitted with a database and a database manager. 
The database is comprised of various measures of Chinese domestic political 
crisis and policy Indicators. Tests of the system's forecast reliability are 
performed with these indicators and the results have verified several 
hypotheses regarding Chinese politics. 

1.1 Research Objectives 

In order to enhance the forecasting capability of the EUAMS system 
this study has undertaken the following objectives: 

1. To develop an Interactive computer forecasting system which, 
with minimum user inputs, will automatically generate in real 
time crisis predictions from a single series of observations or 
through multiple leading indicators. 

2. To construct various Indicators for the purpose of testing 
the forecast reliability of the CEVM system using hypotheses 
linking elite crisis behavior to policy changes In China. 

Attaining the above objectives has enhanced the capabilities of the 
EWAMS system in terms of 

- its computer-based analytic and forecasting capabilities 

- Its use of advanced forecasting methodology 

- its development of crisis Indicators 

1.2 Major Research Accomplishments 

1. Interactive, automated real time computer forecasting system 
A computer forecasting software system has been developed that 
permits direct user-system Interaction, performs automated model 
fitting, and generates Instantaneous forecasts. The system also 


permits the combining and the updating of forecasts. With a minimum 










of user input the system Is largely self-driven. 


2. Indicator Construction 

Various crisis and domestic policy indicators have been 
constructed. These Indicators are developed from elite bio¬ 
graphical data sources and aggregated from economic and defense 
data. 

3. Testing of Hypotheses 

Use of the indicators to test the system's forecasting reliability 
has permitted the verification of two alternative hypotheses per¬ 
taining to factional fractionation in Chinese politics during crisis 
and to the relationship between crisis and policy changes. 



This report is divided into three parts. Part I gives a non-technlcal 
general description of the project. It covers 

- modeling philosophy, relating to design, mathematical background, and 
system testing. 

- computer programming in relation to developing the various functional 
and structural attributes of the system. 

- results of the system’s forecasting reliability tests linking the crisis 
and the policy indicators. 

Part II provides technical summaries of the mathematical background, the 
coding procedure, and the statistical models used in designing and testing the 


CEW1FC system. 






2.0 General System Philosophy 


CEVMFC Is an integrated data management and forecasting system. 

It links a database to a forecaster In an Integrated unit that provides for 
the permanent storage and periodic updating of crisis data for which forecasts 
are generated In real time. For the purpose of ear^y warning and monitoring 
of crises under dynamic conditions, the system is designed to operate in real 
time permitting direct user-system interaction. This interface is facilitated 
by a system of interactive hardware units and software packages. Through this 
conduit of hardware and software the user performs a range of statistical 
and manipulative tasks resulting in short-range multivariate forecasts of 
crisis frequencies and crisis probabilities. 

Moreover. CEUMFC is designed to operate as an automated system that 
requires minimum user operational guidance. All technical modeling decisions 
are pre-progracmed into driver routines directing model fitting and forecasting. 
This feature of the system simplifies the user-system interface and renders It 
accessible to users with limited knowledge of modeling and forecasting metho¬ 
dology. The system so created is thus simple to use, efficient to operate, 
and accessible to technical and non-technical users alike. 

Nonetheless, certain system maintenance must still be performed. The 
inputlng, editing, and updating of data demand some familiarity with programing 
general information systems. It is therefore assumed that the system will 
be supported by a technical staff that maintains the integrity of the database 
and the operational efficiency of the database manager and of the forecaster. 
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2.1 System Design Operational Mode 

The system maintains an array of software packages encompassing over 
100 subroutines, supported by ARPA's Demonstration and Development 
Facilities (DDF) hardware units. The software packages are arranged in a 
hierarchical structure of layers of controlling main programs and computational 
subroutines as Figure 1 illustrates. At the top of this structure is the 
Interface component that elicits, receives, and forwards user commands, and 
that arranges and activates various major components of the system. 

Immediately beneath this level are the software packages subdivided into 
one of two general categories: those that manage the database and those 
that direct the modeling and the forecasting tasks of the system. These 
Include various routines and subroutines that regulate the filing, re¬ 
trieving, securing, and editing of data and execute the model fitting and 
forecasting. 

Program execution is expedited in terms of a branching process commonly 
known as "overlay structure". This process permits the system to focus on 
a sequence of tasks performed by a chain of routines that constitute merely 
a branch in the overall structure such as that enclosed by dotted lines in 
Figure I. By focusing on a branch of the system to complete a particular 
sequence of tasks, this branching process ensures significant saving In space 
allocation which is essential for the efficient use of a minicomputer such as 
the PDP-11/70 of the EVAMS system. 

2.2 The Model 

The modeling and forecasting softwares are built upon a time series 
methodology. A class of discrete-time models known as the "autoregressive- 
moving average vector" (AR.MAV) models form the mathematical foiihdatlon of 
the CEWMFC system. As section A.1 illustrates, these models are systems 
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of extended "autoregressive-moving average” equations describing a structure 
of dynamic interrelationships: 

’t •! A-1 ♦ ?.!V. * *t - !i?Vi - - iA-i 

where the Y t 's are pxl vectors of lagged time series, and a t 's are pxl 
vectors of white noise series, the f's and I's are pxp matrices of auto¬ 
regressive and moving average parameters respectively, and the B d 's are matrices 
of difference operators. These elements are further discussed in section A.l. 

The model assumes the dependence of each series, Y^, on its own pasts, Y^ 
Cor J>0 . and on the past noises, f° r J >0, as well as the dynamic 

Interdependence between series, e.g., Y^'s dependent on Y pt _j* 1/1 other 
words, the model links multiple input and output series via a network of 
causal feedback relationships. The dependence of successive states within 
series represents the memory dynamics that filter through time, whereas the 
Interdependence between series represents their dynamic feedback relationships. 

ttiile memory dynamics from the past are assumed to filter through time 
to the present no reverse process is assumed. For Instance, the model does 


Interdependence between successive noise series. 


2 .3 Modeling Strategy 

The modeling strategy for theCEWMFC system merges a time series modeling 
procedure with a step-vise regression approach. A popular time series modeling 
procedure attributed to Box and Jenkins (1970) prescribes a three-step approach, 
encompassing Identification, estimation, and diagnostic checking. Identification 
is the process by which a tentative model is entertained. Through estimation we derive 
the parameter estimates. Through diagnostic checking we test for model in¬ 
adequacy through analysis of residuals. In contrast, the step-vise regression 






1 
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approach calls for successive approximations of higher order models. We 
thus identify time series models via successive approximations of model to 
data, we estimate models by nonlinear least-squares method, and we check the 
final approximating model by a goodness of fit test. 

The merging of step-wise regression approach to time series analysis 
is essential to the development of an automated modeling procedure, as the 
normal time series identification method Involves rather technical procedures. 
Thus, as in step-wise regression the present approach seeks to explain away 
as much of the variations in time series as possible via successive model 
approximations. Because time series dynamics are caused by temporal de¬ 
pendence between successive observations, temporal dependence can easily be 
successively approximated by higher order models. Since it merely entails 
successively increment In* the model by an order at a time, the procedure is 
mechanical and easily programmable. 

2.3.1 Identification 

Using a recursive search procedure we can design a computer program 
that automatically scans tnrough ascending orders of time series models for an 
optimal model. It starts with the basic first order autoregressive vector process 
and steadily progresses through higher orders of autoregressive-moving average 
vector processes. As the search proceeds through higher order autoregressive- 
moving average vector models by Increments of one in the autoregressive and 
the moving average components, a pair of models is compared at each Increment 
to determine whether the added parameters in the higher order model enhance its 
explanatory power. This relative explanatory power of the models is ascertained 
by comparing their residuals for significant differences. This search procedure 
continues until no difference is observed and the lower order model is retained 
as the optimal case. 


r 





Starting with ARV(l) the model la Incremented to AKMAV(2,1), and then 


to ARMAV(3,2) and so on. Each Increment adds an order to the autoregressive 
and the moving average components. These are known as the "mainline" models, 
because they represent linear extensions from the basic ARV(l) form. There 
are alternative combinations of autoregressive and moving average components 
to each mainline model but they will not be considered until the best fitting 
mainline model is chosen. For Instance,as Figure 3 in section A.2.1 shows, 
alternative combinations to the mainline ARMAV(3,2) model Include ARMAV(3,3), 
ARMAV(4,2), ARMAV(3,1), ARMAV(2,2). After these models are compared and the 
best model chosen, a search for superfluous parameters in the final model 
concludes Che search procedure. Superfluous parameters are identified by 
comparing models in descending orders of moving average parameters, such as 
ARMAV(4,2), ARMAV(4,1), and ARMAV(4,0). 

Since most data are not free of sampling error, a test of significance 
is used in model comparison. The test compares each pair of models to determine 
whether the observed Improvement in goodness of fit between models is 
"statistically significant" or is attributable to random perturbations. The 
test compares the variance-covariance of the residuals for pairs of models. 

The test statistic has an F distribution. As it is shown in section 
A.2.2.4, the residual variance-covariance for our model is best described 
in matrix form (as our model consists of a system of equations which are 
interdependent and autocorrelated). The F test compares the determinants of 
the residual variance-covariance matrices for pairs of models. By setting up 
a critical test value at a given significant level under the null hypothesis 
that there is no significant difference in the models' explanatory powers, 
we can build Into our program a decision rule to compare our F test 
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statistic and to accept or reject a higher order model as a better alternative 
to the lower order case. This search for a better higher order model continues 
until the test statistic reveals a value that, when compared with the critical 
value, suggests acceptance of the null hypothesis and retention of the lower 
order model. Uhen no further improvement is possible the program will auto¬ 
matically determinate all further computation. 

For a fixed critical value it la possible to arrive at several significant 
models, all of which pass the critical level test. But the models may not 
provide identical forecasting results. Thus, an additional test of forecast 
reliability Is also performed for comparison among the alternative significant 
models. This reliability test compares the "mean square error" of the 
forecasts (sum of squares of the differences between the forecasts and actual 
observations of alternative significant models. The best model is one with 
the minimum forecast mean square error. Accordingly, the model so chosen 
not only is statistically significant but also provides the best forecasts. 

Finally, having chosen the most powerful forecast model, we go through 
a final check for superfluous moving average parameters. According to regres¬ 
sion logic, most time series dynamics can be captured by the systematic 
(the autoregressive in the present case) part of a model, and only noise 
variations should exist in the error structure. Thus, it is entirely possible 
to represent a time series by purely autoregressive models or by coupling 
models with high orders of autoregressive parameters to low orders of moving- 
average parameters. By concurrently incrementing the autoregressive and the 
moving average components, it is possible to lead to overfitting later. To 
teat this possibility, a final check for superfluous moving average parameters 
in the selected forecast model is performed. The model is successively 
examined for decreasing orders of moving average parameters. All zero moving 
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average parameters are deleted from the model before forecasting Is attempted. 

2.3.2 Estimation 

A nonlinear least-squares method is used In parameter estimation. Since 
the nonlinear least-squares method uses an iterative search procedure to 
arrive at the final maximum likelihood estimates, it requires the user to 
supply the initial values to start off the iterative search procedure. How 
these initial values are determined affects the quality of the final estimates. 
As it is beyond the means of most general users to supply good initial values, 
we employ a three-step estimation procedure that automatically derives the 
initial values which are ordinarily supplied by the user. The approach 
develops a method whereby the initial parameter estimates for initializing 
nonlinear least-squares estimation is automatically derived within the system. 
The steps are as follows. 

(1) A linear autoregressive vector (ARV(k)) model is first fitted to 
data by means of the linear least-squares method. Estimation is direct 
requiring no user inputs. 

(2) The linear parameters estimated are then transformed, by an inverse 
function, into a corresponding set of nonlinear autoregressive-moving average 
vector (ARMAV(m,n)) parameters. The function links the linear parameters of 
infinite order to linear combinations of nonlinear parameters, as (A.3.2) 
suggests. Equating the linear least-squares estimates to the inverse function 
parameters yields the necessary information for computing the corresponding 
nonlinear autoregressive-moving average estimates. Since the inverse function 
parameters are linear combinations of the autoregressive-moving average 
parameters, it is possible to compute the latter through a recursive procedure 


once the former are known. 
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(3) The transformed autoregressive-moving aversge estimates are then 
taken as the Initial values for starting off the iterations in nonlinear 
least-squares estimation. The nonlinear least-squares estimation method 
uses the Fletcher-Powell (Fletcher and Powell, 1963) minimization method to 
minimize the realdual sum-of-squares matrix of the vector model. More detailed 
technical discussion is given in section A.2.4. 

2.3.3 Diagnostic Checking 

A check for goodness of fit of the selected forecast model concludes the 
modeling procedure. This final stage further tests the adequacy of the model. 
It Involves a test for systematic trends in the residuals of the fitted model. 
The presence of such trends indicates model inadequacy. Several statistics 
are employed at this stage. They include the residual autocorrelations and 
partial autocorrelations, the cross-correlations between the residuals and 
the time series, and the associated Q and S significance tests. These 
statistics address not only trends within series but also trends between 
series. 

2.4 Forecasting 

The objective of time series modeling is the development of optimal 
forecasts.These forecasts produce th4 least amount of forecast mean square 
error in relation to all possible alternative forecasts. 

These minimum mean square error forecasts are recursively computed 
from the forecast equations according to some rules of conditional expec¬ 
tation. The results contain both elements of naive autoregressive-moving 
average forecasts and leading indicator forecasts. The CEVMFC system thus 
can generate both naive forecasts from a single series or forecasts by 
leading Indicators with several series. 
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3.0 The Computer Forecasting System 

The CEWMFC system is programmed in ANSI standard Fortran IV language. 
The system operates In real time to permit direct system-user Interface. 

Tt> simplify system operation, all necessary user command Inputs are prompted 
by system solicitations, which appear in the form of questions, such as 
"number of forecasts?" This streamlines user Input commands in relation 
to a prearranged sequence of operations. Moreover, the user Is also aided 
by "help" routines which provide additional explanations of and instructions 
for the actions required. As moat technical decisions are Internalized In 
the driver routines, the solicitations are limited to simple, nontechnical 
decisions. 

3.1 System Functional Attributes 

As a management aid for crisis monitoring, the CEWMFC system generates 
crisis frequencies and crisis probabilities. Its forecasting capability Is 
enhanced by the following system attributes. 

3.1.1 Interactive and Real Time 

The system Is designed to Interface directly with the user via graphics 
terminals. User commands induce lmedlate responses from the system. This 
Is particularly Important for solving sequential decision problems that 
frequently confront the crisis analyst. As crisis events unfold In rapid 
succession and as new data become available with each unfolding event. It 
Is Imperative that forecasts be generated quickly In order to ensure 
continuous monitoring of events. Thus, the ability of the system to Interact 
In real time with the user Is an essential feature of the CEWMFC system. 

The automated modeling capability also enhances the system and Its forecast 
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updating capability which further minimizes user decision time and ensures 
rapid fine tuning of previous forecasts. 

3.1.2 Automated Modeling 

While most existing time series techniques (e.g. the Box-Jenkin's (1970) 
approach or Brown's (1966) exponential smoothing) require the user to Identify 
the correct time series model, the CEVMFC system uses an automated procedure 
to identify the optimal forecasting model. Through a combination of significance 
test and forecast mean square error test nested in a recursive model search 
routine, the CEUMFC system automatically identifies and estimates the optimal 
model. Such technical decisions normally required in model building are in¬ 
corporated in driver routines. Thus, relieved of the tasks of making technical 
decisions, most non-technically trained users may find the system simple to 
operate. 

In addition, some optional statistical maneuvers are available for the 
trained specialists. These statistics permit the user to manipulate complex 
time series (e.g., nonstationary series) using various transformational 
procedures and to check for the adequacy of the optimal model determined by 
the system. 

3.1.3 Multiple Forecasting Capabilities 

CEWMFC system maintains multiple forecasting capabilities. The 
system generates both "projective" and "objective" forecasts (Andriole and 
Young, 1977). Projective forecasts are trends extrapolated from single series 
of observations. Objective forecasts are forecasts generated by leading in¬ 
dicators. These forecasts may be used interchangeably to address a variety 
of research needs. Besides providing naive numerical projection, they may be 
used to study Interaction pattern*., systems and stability analyses, and 
optimal control analysis. Some of these potential applications are discussed 






15 


in a series of monographs by Zinnes, Gillespie and their co-workers 
(1978a, 1978b. 1979). 

The system updates forecasts with each new piece of Information using 
a simple "adaptive forecasting" method which does not require refitting 
the entire forecasting function. The method adjusts for the forecast error, 
revealed by the latest piece of Information entering the system, through 
a simple weighting function (Section A.3.4). The result is a continuous 
series of adjustments that fine tunes the forecasts over time. Such 
adjustments are localized about single pieces of Information; they are not 
sensitive to abrupt changes or to steady long-term drifts In the time series. 
Updating must therefore be limited to short-range forecasting. 

The system maintains combining forecasts capability. Forecasts gen¬ 
erated by several alternative methods are combined when they are available 
for a given Indicator. The combined forecasts have been shown to yield more 
accurate results than their original values (Dickinson, 1975). The method 
Incorporated In the CEVIMFC system combines several forecasts, each weighted 
by its estimated sum-of-squares of past forecast errors. The forecasts to 
be combined are not restricted to those generated within the system. Fore¬ 
casts generated externally may also be Incorporated Into the system. The 
system thus serves an integrative function In linking forecasts from alternative 
forecasting systems. 

For each frequency forecast the system computes its associated 
probability. The probability Is derived from a normal distribution and it 
measures the likelihood of the forecast to deviate from an average of past 
values by Its observed magnitude if only random effects were present. Abnormal 
deviations Imply crises. 
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3.1.A File Manipulation Capability 

The CEVMFC system maintains a database manager which ensures the 
accessibility, integrity and security of data entered into the system and 
stores the data generated within the system. The system receives incoming 
data from external sources; it creates new data using various transformations; 
it edits existing data in accordance with user commands. 

The system maintains four separate file types, including the directory file 
tike descriptor files, the data files, and the model files. These files 
identify the countries for which data are stored, store the data, describe 
the content and arrangement of the data stored, and store the statistical 
results of model fitting and forecasting. 

Moreover, the system is capable of storing and manipulating data for a 
maximum of 100 countries and for each a total of 100 indicators. For each 
indicator the system may maintain a maximum of 500 observations. Finally, 
the system has the capability to accomodate data for a total of five estimated 
models. 

Most types of files, excepting the data files, are programmed in 
unformatted direct access (random) form that permits direct user access to 
any part of the file structure, eliminating the need to search through a 
prearranged sequence of information as in formatted sequential files. 

Such a flexible format structure shortens operational time required in 
locating data within the system. Nonetheless, no direct access to the files 
independent of the system is available to the user. Since the system 
requires a sequence of executions to effect data editing and daca deletion, 
it does secure the data against unintended alteration or destruction. 

3.2 System Structural Attributes 

The CEVMFC system is comprised of three main structural components , 
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Figure 2 

CEWMFC System Structure 
















as Figure 2 illustrates. They are 


a. A set of command modules which monitor the lnputlng, the outputlng, 
and the editing of data; the modules regulate the graphic displays 
of the forecasts; and control the implementation of statistical 
analysis to identify, estimate, and generate forecasts 

forecast functions. These modules are tied together by a main 
program which is the interface that facilitates the real time 
interaction between the system and the user. 

b. A database manager which regulates the filing, retrieval, securing, 
editing, and storage of data. 

c. A system of statistical routines that Identify, estimate, and generate 
forecasts from the forecast functions as well as perform various 
auxilllary statistical computations. 

3.2.1 The Control Modules 

The control Modules consist of four components tied together by a main 
program that serves as the Interface between the system and the user. This 
interfacing routine receives and deciphers incoming user commands and 
activates various components of the system. Subsumed under this main program 
are four control modules: 

(1) A data entry control module that regulates the lnputlng and the 
outputlng of data. This module controls the creation of data 
and its entry into the system. Data creation is achieved through 
one of the ensuing transformations: 

-logarithmic 

-multiplicative, indicator by scalar 
-multiplicative, indicator by indicator 


-exponential power 
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-linear combination 
-regular difference 
-seasonal difference 

(2) A data editor/lister that lists data and system information 
(e.g., name of country and its associated descriptor file) 
and edits such information as the user specifies. 

(3) A data display control module that regulates the graphical and 
tabular outputs from the system. 

(4) An analytical control module that coordinates the executions of 
an Integrated system of time series model fitting and forecasting 
algorithms. 

3.2.2 Database Manager 

The database manager includes a set of programs that maintains the 
security and organizes the storage of data. The data Include raw and trans¬ 
formed observations, the forecasts and their associated statistics, and the 
parameter estimates. This manager monitors user access to the data file, 
maintains the Internal consistency of data, secures the data against lnadver- 
tant alternation or destruction, and organizes the filing and storage of 
data. The databse Itself is subdivided Into a core and a file structure. 

3.2.3 Analytic Routines 

The analytic routines form an integrated system of algorithms which 
Identify, estimate, and check the goodness of fit of the forecast functions, 
and generate optimal forecast as well as various associated statistics. These 
algorithms are arranged hierarchically to perform the following sequence of 
operations. 


(t) Optional autocorrelation, partial autocorrelation, and cross¬ 
correlation functional analvsls permit the user to scan for evidence of 
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nonstationarlty In the time aeries and to Identify dynaalc 
relationships between series. These statistics provide reasonable 
clues to time series nonstatlonarity, as evidenced by non-declining 
autocorrelation and partial autocorrelation values (Box and Jenkins, 
1970). But only a rough approximation of the underlying dynamic 
relationships between series can be expected from cross-correlational 
analysis, as cross correlations are inefficient estimates con¬ 
founded by autocorrelation effects within series. The Nonstationary 
time series may be rendered stationary by a difference or a 
logarithmic transformation. These transformations are effected by 
calling the data entry module. 

(2) Assuming stationerlty, the parameters of an autoregressive vector 
model are first estimated by linear least-squares method. 

(3) The autoregressive vector estimates are then transformed into 
nonlinear autoregressive-moving average vector parameter estimates 
via an Inverse function. The Inverse function links the linear 
autoregressive parameters to linear combinations of nonlinear 
autoregressive-moving average parameters. 

(4) These transformed nonlinear estimates are then entered as initial 
values into a nonlinear least-squares estimation routine. These 
values start off the Iterations in the nonlinear estimation method 
which eventually leads to the maximum likelihood estimates. 

(3) Next, autocorrelations and partial autocorrelations of the residual 
are derived for the fitted model, and cross-correlations between 
pairs of residual series and actual observations are computed. 

These statistics should reveal any remaining stochastic trends 
within each time series or any causal relationships left unexplained 
by the model. In addition, chi-square tests of 






significance are performed to differentiate significant 
trends from random fluctuations. When no significant 
trends are revealed by the autocorrelations, by the partial 
autocorrelations, and by the cross-correlations, the fitted 
model is adequate. 

(6) Parameter estimates from the fitted model are transmitted to 
a main forecasting program which, in turn, determines the form 
of forecasts to be generated. This main program coordinates 
several forecasting algorithms. 

(a) A general algorithm which generates multiple series of 
forecasts on the basis of past autoregressive-moving 
average dynamics and of the leading indicators. Moreover, 
the algorithm adjusts for all difference operations performed 
on the original series by re-incorporating them into the 
forecasts. 

(b) An algorithm computes the variance of the forecasts and 
constructs the confidence bands for each forecast. 

(c) An adaptive forecasting algorithm updates existing fore¬ 
casts by adjusting for the forecast error revealed by each new 
piece of information entering the system. 

(d) An algorithm that combines forecasts generated by 
alternative methods to arrive at a single series of forecasts. 

(e) An algorithm computes the standard normal probability 
associated with each forecast. 

4.0 System Testing 

Tests of validity and forecast reliability are Implemented at each step 
of CEUMFC system construction. Each subroutine is designed to achieve a 
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certain performance level and la tested at each step of system development. 

The test employed is known as "sensitivity analysis" (De Neufvllle and 
Stafford, 1971). It gauges how a change in one part of the system effects 
the whole. By this method either we vary the parameter values or we vary 
the input series and Inspect the corresponding changes in the output series for 
evidence of possible instability. 

Moreover, upon completion of the entire system, a test of forecast 
reliability is performed with data secured within the system. 

4.1 Forecast Reliability Test 

The test, known as "retrospective forecasting", gauges the historical 
accuracy of the forecasts. By this method, a segment of the latest data 
is compared with forecasts for the same time period. Large deviations 
of the forecasts from the actual observations would cast doubt on the 
reliability of the forecast system. 

Since random disturbances are always present in the forecasts, dis¬ 
crepancies are inevitable. In order to differentiate systematic from 
randomly Induced deviations of forecasts from actual observations, a 
significance test is Introduced. The test Statistics, devised by G.E.P. 

Box and G. Tlao (1977) is: Q - e'V *e which is distributed as chi squares with 
a degrees of freedom, and where e is a 1 x a vector of forecast errors de¬ 
fined as the difference between the actual observations and the forecasts, i.e., 
e • Z - Z(l) and V 1 • T'T/o 2 , where? is an a x m aatrlx of coefficients defining 
the relationships between time series white noise processes and forecast errors 
generated from the same process, e.g., a £ - Te^ The test addresses the null 


hypothesis that the observed deviations are results of random disturbances 
and not indicative of systematic errors. 
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This test has a major drawback, for It uses the same data that guided 
the formulation of the model to test forecast accuracy. The result is a 
tendency to understate the amount of forecast errors and thus bias the outcome 
in favor of confirming forecast reliability. 

4.2 Hypotheses 

For testing purpose, indicators are constructed specifically for the 
CEUHFC system. These indicators include various political crisis and 
economic/defense policy measures. The selection criteria draw heavily upon 
current social science theories which may be stated in hypothetical form. 

Two hypotheses underlie the selection of indicators. These hypotheses 
adopt the factional model in elite analysis. By definition elites manipulate 
power, and power or inter-factional power distribution is the focal point of 
our model. The model analyzes factional groupings to determine how behavioral 
characteristics and career mobility patterns of faction members are associated 
with political crises and policy outcomes. Briefly, these hypotheses may be 
stated as follows. 

H^: The frequency of purges and demotions Increases during periods 

of political crises and displays patterns of collective ex¬ 
perience shared by members of elite factions set apart by 
competing Interests and value orientations. 

H,: Changes in leadership resulting from factional differences 

lead to significant readjustments in policy priorities. 

Underlying these hypotheses are the assumptions that politics In China 
since 1949 has been largely dominated by factional power struggles and that 
while policy differences do play a role they serve more often as the pretexts 
for rivalry than as the causes of conflict. Nevertheless, once a conflict is 
resolved the dominant faction re-determlnes the policy priorities. The logic 
behind these assumptions is developed as follows. 

With "politics in command" (Gray, 1974) of all policy decisions regarding 
production, the distribution of wealth and authoritative allocation of values. 







24 

differences Invariably will emerge to divide the elites along diverging lines 
of interests and value orientations. Along these cleavages, elite factions 
tend to coalesce about natural alliances based on shared background attri¬ 
butes and common career experiences (Chang, 1969). From such collective 
experiences, there emerge the personal ties essential for the working 
of factional politics. Nathan (1973) has observed that the factions, principally 
at Che central governmental level, form the bases through which personal 
ambitions are fostered and political power accrued. As the power balance 
tips in favor of partioular factions, groups of elites share collective 
experiences of career successes and failures. Thus, a recent Rand study 
(Sung, 1975) found that during crisis periods, elites with certain shared 
interests and background attributes are more likely to have similar occu¬ 
pational mobility patterns than those who do not. 

Since elite factions are defined along the lines of existing power 
cleavages, the historical twists and turns of their relative power positions 
offer illuminating insights into trends of power distribution in Chinese 
politics. As the power position of each faction depends largely on the 
extent of control its members can exercise on Important party and governmental 
decision-making organs, the rate of positional loss such as demotion has a 
profound effect on prevailing power distribution. While a normal rate of 
personnel loss resulting from natural causes is Inevitable, a high turnover 
rate in an otherwise stable elite group is indicative of changing factional 
power balance. Accordingly. Oksenberg (1974) has found in his analysis of 
the Chinese elite "exit pattern" that a high rate of purges and demotions is 
the usual pattern of personnel attrition in times of crises. As major 
power redistribution among factions invariably disturbs the leadership It 
creates stresses in the political system and threatens to alter the existing 
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political objectives and policy priorities. The resultant uncertainty is 
likely to cause a cessation of normal activities in the political process. 
Thus, while the relative personnel loss rates among factions are not direct 
measures of political crises, they are nonethless good proxy indicators. 

Since power redistribution among factions in the leadership often 
leads to policy changes. It carries profound policy implications for Chinese 
politics. Most great movements of the recent past, the Crest Leap Forward 
of 1958-1960, the Cultural Revolution of 1966-1969, and the purge of the 
"radical four" in 1976, have been marked by high levels of personnel losses, 
followed closely by major policy changes. Tracing the broad configurations 
of historical evidence, scholars have fashioned various explanations of 
Chinese politics. While all have observed some oscillations in patterns of 
policy changes some have attributed those patterns to impacts of crises re¬ 
sulting from "left" versus "right" and "idealogue" versus "pragmatist" power 
struggles (Nathan, 1976). Parris Chang (1975), for example, notes that policy 
shifts in China arc often associated with changes in the power balance 
among contending "left-right" factions in the decision-making councils. 

The linkage between political crises and policy changes is supported 
by a rich heritage of historical and empirical evidence (Rosenau, 1969; 

Liao, 1976; Nathan, 1976). In domestic policies, one may trace a clear 
chronology of parallel movements in crisis developments and policy shifts. 

For example, the Creat Leap Forward of 1958-1960 initiated an era of mass- 
mobilisation culminating in policies of collectivization and rapid In¬ 
dustrialization. The excesses of the Great Leap then led to a period of 
retrenchment towards decentralization and tolerance of material Incentives 
lasting from the late 1950's to i965. With the coming of the Cultural 
Revolution in 1965, all pragmatic policies associated with planning. 
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managerialism, specialization and technology were cast to the wind, to 
be superceded by ideological dogmatism. After four tumultuous years, 

Chinese policies since 1969 have regained some degree of moderation — in 
spite of a brief Interlude of radicalism in the mid-1970's — and redound 
to a search for pragmatic solution. These events, among other historical and 
empirical evidence, supply substance to a linkage between crises and policy 
changes in China. 

In short, it is suggested by the above discussion that (a) factions 
embody conflicting interests and value orientations and the relative rates of 
personnel losses between factions provide good indications of political crises, 
and (b) political crises Induce changes in Chinese military and economic 
policies. 

4.3 Data 

Data are collected for various Chinese political crisis and economic/ 
military policy indicators. For the crisis indicator a measure of elite 
turnover rates is used. The turnover rates are measured in terms of the 
frequency of purges and demotions. 

Although there is no totally satisfactory way of identifying factions, 
one possibility is to define factions along biographical attributes that 
imply shared values and Interests. For example, the following attributes 
are selected for this study. 

-commander versus conlssar career patterns 
-military versus civilian career patterns 
-administrative "functional" affiliations 
-military region and field army affiliations 
-generational affiliations 
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Aa defined along these biographical attributes the relative turnover 
rates within the factions reveal the pattern of inter-factional conflict 
signaling the onset of crisis. 

Approxlnately 568 top level party and governmental personnel are 
selected for the analysis, including members of the Politburo and of the 
Central CoMittee, ministers of the state Council, commanders and first 
political commissars of the military regions and the service arms, and the 
governors and first party secretaries in the provinces and on the revolutionary 
committees. Of these, certain biographical data are available for 400 elites 
as compiled in a previous Rand elite biographical study covering the period 
1956 to 1973 (Sung, 1975).* Our analysis thus extends the data collection 
effort through 1978. 

Since the Central Conmilttee has undergone substantial changes in 
membership since 1974, 86 new members have been Identified. For these 
elites we have also collected various other biographical data which are not 
used in testing the present system but which are nonetheless necessary for 
maintaining some consistency in extending the earlier data set. These 
attributes obtained for the new members are: 

-place and date of birth, 

-education, military and civilian, 

-combat experience, 

-awards, 

-date entered party. 

Annual data for the policy indicators are collected along different 
economic and military dimensions from 1956 to 1978. They Include: 

-gross national product, 

-export and import levels, 

*Two cases, Yang Ching-Jen and Chlang Wel-chlng, were each Inappropriately 
coded twice aa separate individuals in the Rand study. Correcting for 
these errors reduces our sample size by two. 
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-total industrial outputs, 

-total agricultural outputs, 

-defense expenditures. 

It should be noted that the quality of Chinese statistics is limited 
at best. For a combination of reasons, including a traditional lack of 
concern for numerical accuracy, the manipulation of data for reasons of 
political expediency, the infancy of the Chinese State Statistical Bureau, 
the central data collection agency, and the disruptions of the Creat Leap 
forward and the Cultural Revolution to statistical research efforts in 
China, statistics for most years are not totally accurate and exist for 
some years only in the form of rough estimates (Orleans, 1974). Nonetheless, 
the aforementioned data are adequate for the limited present purpose of 
testing the forecast reliability of the CEWMFC system. 

4.4 Indicator Construction 

Two sets of indicators are constructed to drive the CEWMFC system. 

4.4. I Crisis Indicators 

Several steps are required to construct the crisis indicator. 

First, we assign our sample of elites to various factions sharing common 
biographical attributes. Thus we have identified various factions according 
to their party affiliations, their governmental functional affiliations, 
their field army affiliations, and their generational affiliations. 

Second, we compile the total annual turnover rates due to purges or 
demotion for the factions. 

Third, in order to capture the dynamics of inter-factional conflict, 
the turnover rates are further transformed into ratio scores. These scores 
are obtained by Juxtaposing competing factions and dividing the turnover 
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frequency of one Into that of the other. For Instance, 

number of purged and denoted commanders 
i number of purged and demoted commissars 


Y, 


number of purged or demoted personnel In field army i 
number of purged or demoted personnel In field army J 


m number of purged or demoted personnel In administrative functional area 1 

j number of purged or demoted personnel In administrative functional area J 

number of purged or demoted personnel in generation 1 

4 " number of purged or demoted personnel in generation j 

Y ■ number of purged and demoted civilians 

number of purged and denoted military personnel 

Finally, because some personnel loss will always occur in any given 
faction from natural attrition, it is necessary to differentiate crisis rate 
from the normal loss rate. An index of abnormality, the Z score, is used for 
this purpose and it is measured according to 



where Y f is the average score for past Y t 's and the corresponding standard 
deviation. As Andriole and Young (1977) noted, this score is akin to the 
standard normal score and it measures the degree of departure of a current 
value from the past average. A large Z score suggests abnormal departure from 
a state of "usualness" as Che average is caken to imply. 


4.4.2 Policy Indicators 

The policy indicators are similarly computed in terms of ratio scores 
to tap Che conflict dimension involving competing Issues in the policy-making 
process. Again, Juxtaposing measures of competing policies, we take the 
ratio of their respective scores to produce the following indicators. 
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^ m industrial outputs 
1 agricultural outputs 

^ _ defense expenditures 

T 2 gross national product 

^ m total import _ 

3 gross national product 

y m total export _ 

A gross national product 

These ratios are Intended to reveal the alternating swings between competing 
policy positions taken by the Chinese leadership between 1956 and 1977. 

A•5 System Testing 

The CEWMFC system Is tested along three dimensions. Sensitivity analysis 
assesses the stability and the responsiveness of the system to parameter 
changes. Input variations, and effect of perturbations in subcomponents on the 
overall system. Forecast reliability test ascertains the accuracy of the 
forecasts generated by the system. Finally, hypothesis testing tests the 
ability of the system to discern complex network of Interrelationships to 
generate the corresponding parameter estimates, and to provide multiple 
forecasts. 

A.5.1 Sensitivity Analysis 

Sensitivity analyses at every level of the CESMFC system reveal the 
system on the whole to be stable. With further explication provided by 
"help" statements, the system prompts user consnands with direct requests 
for specific Inputs. The choices of Inputs requested are displayed on 
screen, thus restricting the user to a limited range of operations. All 
unacceptable commands are screened out by the Interface routines, and the 
correct commands are suggested by help statements and error messages. This 
has the effect of Insulating the system from Inappropriate command Inputs. 
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The database manager Is, in contrast, totally insulated from external 
inputs. No direct user access is permitted into the manager, and since 
the manager also operates independently of data, it is a completely self 
sufficient unit. As a result, this component of the system is highly 
stable and functions according to a routine chain of operations. The 
computational algorithms remain largely stable over a wide range of 
parameter values. Parameter changes in most subroutines do not adversely 
affect the performance of the algorithms. This holds crue for the linear 
least-squares routine used in estimating the autoregressive vector model, 
the inverse function routine, and the various forecasting routines. The 
exception is the nonlinear least-squares routine for autoregressive-moving 
average vector model estimation. Unlike other routines, the routine is 
sensitive to parameter perturbations. Although a superior modification of 
the gradient search method, the Fletcher-Povell optimization algorithm used 
in the nonlinear least-squares routine is limited in the number of parameters 
which it can effectively estimate. While the Fletcher-Powell method is to some 
extent data dependent, in that model overfitting impairs its estimation 
efficiency, it has been shown to be capable of accommodating at most five 
series at a time in lower order vector models. As the order of the model 
successively Increments, the number of iterations required for convergence 
correspondingly Increases. Our experience has shown that in particular for 
overfitted model with large number of parameters, convergence is difficult to 
achieve Many iterations required in computing the optimum step size used 
in determining the direction vector and in achieving convergence towards the 
maximum likelihood estimates. < Introducing some termination rules, stopping 
the iterations after a certain number has been reached, has served to hasten 
convergence. But the effect on the estlout ion may be negative especially 




In situations where the likelihood function is not quadratic and well behaved. 
Table 1 shows the number of Iterations required before convergence convergence 
is achieved for different orders of models, using one and two series of 
observations. 

Moreover, sensitivity analyses have revealed that due to 
the limiting numerical capacity (small word size) of the PDP 11/70 minicomputer, 
the CEWKFC system can not efficiently estimate models with a large number of 
parameters. Testing with five series at a time often results in floating 
overflows at several points of the system. The determinants of matrices, such 
as che residual sum-of-squarcs-and-cross-product matrices, for higher order 
ARMAV(m,n) models often assume sizes beyond the limits of the PDP memory 
word capacity (upper limit of 1.7 x 10^®) causing the system to abort 
estimation. Writing the program in double precision form has ameliorated 
Che problem but has also doubled the space requirement of the system in 
relation to single precision format. Furthermore, setting upper bounds to 
limit determinant size has also been attempted. Such bounds prevent the 
system from aborting in extreme cases of large determinant values but they 
also impair the system's ability to generate maximum likelihood estimates. 
Nonetheless, since large determinants of residual matrices often are 
indicative of ill fitted models, which should be eliminated from consideration 
by the model search routine, generating inefficient estimates for those 
models should not affect the selection and estimation of the final optimal 
model. The optimal model will in all likelihood produce small residual values. 

Moreover, use of difference operation to remove nonstationarity in time 
series is highly recommended. Estimation with nonstationarity series seldom 
leads to rapid convergence and less frequently to efficient estimates. When 
the nonstationary parameters have values close to 1 differencing has the 
added advantage of reducing the number of parameters that have to be estimated. 
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DUMBER OF ITERATIONS FOR CONVERGENCE 
OF DIFFERENT ORDERS OF ARMAV(a.n) MODELS 


Series 


Orders 


( 1 . 0 ) ( 2 . 0 ) 


(1.1) (2.1) (3.2) (4,3) 


Chinese 
cosaander/ 
coonissar 


So. of 

Iterations 1 27 50 73 


F Test 

significance 

level .05 .05 .05 .05 


Gertsen 

Saval 

Expenditure 


So. of 

iterations 1 22 51 35 


F Test 

significance 

level .05 — .05 


Iterative 

Anglo-O.eraan 

Saval 

Expenditures 


So. of 

iterations 41 50 51 82 


F Test 

significance 

level .05 
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This would hasten convergence and lessen the likelihood of overflowing system 
capacity. Testing using different data sets. Including the Chinese crisis 
and policy Indicator series, monthly Wisconsin employment series, annual 
Anglo-Cerman naval expenditure series (1870-1914) have revealed the system 
to be particularly sensitive to changes In nonstationary series. It Is thus 
recommended that all series entered Into the system be stationary or be 
rendered stationary by means of the difference or the logarithmic transformation 
In the system. 

Similarly, the system any be sensitive to changes In time series when 
the number of observations Is small. Since In short series, each observation 
can be expected to contribute more to determine the nature of the data 
distribution than would an observation in a lengthier series, minor variations 
In a few observations can be expected to affect and alter the forecast outcomes. 
Increasing the number of observations In time series should help to lessen 
system sensitivity to minor variations In time series. 

Finally, the svstem appears to operate best with lower order models In 
which the number of parameters Is limited. Tests of Chinese elite data reveal 
that the system functions most efficiently with less than five series. Con¬ 
vergence towards maximum likelihood values Is rapid and forecasts most preclae 
with single series of observations. The system operates with relative efficiency 
for two and three scries, but encounters more difficulty as It advances towards 
estimating five or more series. Because of the large number of parameters 
Involved In multivariate time series models each successively higher order 
Increases the number of parameters substantially, (e.g., for models with five 
series increment of each order adds 25 parameters to be estimated: ARV(l) 
contains 25 parameters, ARMAV(l.l) contains 50 parameters). Consequently, 
as the number of time series Increases the order of models which the CEWMFC 
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system can effectively estimate becomes progressively smaller. Tables 2a to 
2d show the Inverse relationship between the orders of models and the number 
of series. 

Tables 2a through 2d also Illustrate the automated sequence of 
operations performed by the CEUKFC system as it attempts model fitting. 

Table 2b. for example, shows that the system first examines and compares 
progressively higher orders of mainline models from which It Identifies 
ARMAV(2,1) to give the optimal fit. Then, comparing this model to alternative 
intermediary models. It found only ARMAV(3,l) model to yield significant 
Improvement In relation to ARMAV(2,1). Since, there are several 
statistically significant models, the system attempts a further forecast 
mean-square-error test to single out the best forecast model. The results 
favor retention of ARMAV(3,l) as the final forecast model. 

4.5.2 Forecast Reliability Test 

The forecast accuracy of the CEVJMFC system Is assessed by means of a 
retrospective forecasting method. The method, as described In section 4.1, 
evaluates the performance of the system by the amount of forecast errors 
Incurred. The errors are obtained by comparing retrospective forecasts 
for past periods with corresponding known observations. The forecast errors 
are then examined for systematic, nonrandom discrepancies. Large discrepancies 
in forecast errors are indicative of a lack of reliability In the forecast 
results. 

Results of retrospective forecasting and the attendant statistical 
significance tests for two separate time series are displayed In Tables 
3a and 3b. The columns of the tables display the actual observations, the 
1-step-ahead forecasts from a fixed time origin, the associated forecast 
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errors, Che one-step-ahead forecasts generated from successive time origins, 
the noise ® c '» which are equivalent to the one-step-ahead forecast errors, 
and the squared a^s. The displays also Include the variances of a t and 
the chi-square Q statistics used to determine the statistical significance of 
the forecast errors. 

The test results Indicate that the forecasts are generally in 
agreement with the actual observations in terms of magnitude and of direction 
of movement over time. Creater accuracy seems to characterize the forecasts 
from the employment series than from the Chinese Industry/agriculture indicator 
series. The former generates smaller forecast errors and its movement more 
closely monitors and reflects the movement of the actual observations, 
whereas the latter follows less closely the movement of the actual observations. 

A comparison of the forecast errors reveals that while no distinguishable pattern 
of errors is discernible in the employment forecasts, forecast errors for the 
Chinese series do get larger over time. These divergent results arise in 
part from the fact than the employment series has five times as many 
observations as the Chinese indicator series, which would have contributed 
to more efficient modeling and estimation and thus to more accurate forecasts 


in the former. 
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Table 2a 

Automated Model Fitting Outpu t 
for One Series 


Y £ • Industry/agriculture outputs 


Model 


AHMA(l.O) 
AKMA(2,1) 
ARMA(l.l) 


significant 
not significant 
significant 


location 


mainline 

mainline 

intermediary 


Forecast mean-square error test reveals best model to be 
Ai01A(l,l) 

Final model - ARMA(1,1) 

AR(l) - .62679 
MA(1) —.71071 
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Table 2b 

Automated Model Fitting Outputs 
for Two Series 


Y u - import/CNP 

Y 2t “ generation 1/generation 2 


Model 


F Teat 


location 

ARMAV(1,0 




malnline 

ARMAV(2,1) 


significant 


mainline 

ARMAV(3,2) 


not significant 

if 

mainline 

ARMAV(2,2) 


not significant 

if 

Intermediary 

ARMAV(3.1) 


significant 


intermediary 

ARMAV(2.0) 


not significant 

if 

intermediary 

ARMAV(1. 1) 


not significant 

if 

Intermediary 

mean- 

-square-error 

reveal best model 

Is ARMAV(3, 

,1) 

ARMAV(3,0) 


not significant 

if 

test for zero 





average 

Final model - 

ARMAV(3.1) 




2.0232 

-.29534~ 

p 

ARV(l) 

1.033 

-1.1852 

*_ .51865 

.0312251 

J 

1.2205 

.5964 2 

-1.8559 

-.078638' 




.-.61496 

.45622_ 

ARV(2) 



*.71867 

-.9498f 

ARV(3) 



.. 52372 

. 30357J 





moving 


MA(1) 
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ARMAV(l.O) 
ARMAV(2,1) 
ARMAV(3,2) 

ARMAV(2,2) 
ARMAV(3.1) 
ARMAV(1,1) 
ARMAV(2,0) 


significant 
significant 
not significant 

not significant 
not significant 
not significant 
not significant 


Intermediary 
Intermediary 
Intermediary 
Intermediary and 
lower order moving 
average 


Comparison of forecast mean-square-errors of significant models reveal 
ARMAV (2,1) to yield best forecast. 

Final model - ARMAV(2,1) 

(The CEWMFC system then employs the final model for forecasting purpose). 



Y lt-1 

Y 21—1 

Y 3t-1 

Y lt 

.65716 

-.36788 

. 21350~ 

Y 2t 

-.66100 

-. 10157 ' 

3.7925 

Y 3t 

-1.2828 

.00938 

2.6144 


Y lt-2 

Y 2t-2 

Y 3t-2 

Yu i 

. 58059 

-.35158 

.40317’ 


! -.26823 

. 11792 

-2.3294 

Y 3t 

-.63255 

.14310 

.80377_ 


a 

lt-1 

a 

25-1 

a 

3t-l 

V 

|". 68273 

-.69998 

.74970 

Y 2t 

.85698 

.59884 

.20334 


3t 


2.0397 


.80272 


.40041 
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Table 2d 

Automated Model Fitting Output 
for Four Series 

Y 

It • Industry/agriculture 
Y 2t - lmport/GNP 
Y 3t - export/GNP 
Y 4t - civillan/milltary 

Model F Teat 


ARMAV(1 

.0) 


significant 


ARMAV(2 

.1) 


not significant 


ARMAV(l 

.1) 


not significant 


ARMAV(2 

.0) 


significant 


Forecast mean'square 

test reveals 

best model to be ARMAV(1,0) 

Final Model - ARMAV(1 

.0) 




Y 1C-1 

Y 2t-l 

Y 31—1 

Y 4t-1 

T “ ! 

^"-90788 

.26891 

-.03007 

.048038 

V 

2.4836 

. 27589 

.51484 

.45307 

T n 

1.4570 

.57049 

.76769 

.32054 

T »t 

-1. 7077 

-.0377 

.19815 

. 13675 


location 

malnllne 
mainline 
Intermediary 
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Table 3a 


Forecast Reliability Teac 
Wisconsin Employment Series 
n-104 


steps 

ahead 

Wisconsin 

observations 

series n*104 

i-step-ahead 

forecasts 

1 

550.0 

548.35 

2 

544.9 

538.15 

3 

542.4 

533.45 

4 

539.0 

533.88 

5 

532.5 

528.46 

6 

532.5 

528.46 

7 

533.5 

531.37 


forecast 

one-step-ahead 

a t 

2 

a 

errors 

forecasts 


t 


1.7 

549.15 

.849 

.721 

6.8 

539.87 

5.026 

25.262 

9.0 

541.22 

1.175 

1.380 

5.1 

541.73 

-2.735 

7.485 

4.1 

529.78 

2.512 

6.314 

4.8 

533.52 

-.323 

.104 

2.2 

536.07 

-2.578 

6.648 


Sum 


3.926 47.914 


S 2 - [Eat 2 /n] - C(Ia t ) 2 /n 2 ] - 6.53 
Q - 1/S 2 Ia t - 1/6.53 47.91 - 7.9 

Since chi-square critical level is X 2 ^ - 14.07. test suggests no systematic erroi 


* the a t ’s are the one-step ahead forecast errors computed from 


where Y t (l) is the one-step-ahead forecast from time origin t corresponding 


t+I. 


to the actual observation Y 
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Table 3b 


Forecast Reliability Test 
Chinese Industry/Agriculture series 
n-13 


step 

ahead 

observations 

1-step-ahead 

forecasts 

forecast 

errors 

one-step-ahead 

forecasts 

!% ‘1 

1 

3.3899 

3.1965 

.1934 

3.196 

.193 .0372 

2 

3.9199 

3.0269 

.8830 

2.931 

.988 .9761 

3 

4.2799 

3.0953 

1.1846 

2.973 

1.306 1.7051 


Sum 


2.467 2.7184 


S 2 - [2* t 2 /nj - [(Ia t ) 2 /n 2 ] - .2299 

Q - l/S 2 I « t 2 - 1/.2299 x 2.7184 - 11.787 

Since chi-square critical level Is X 2 q, ^ • 7.81, test suggest presence of 
systematic error. 

•the a c ’* are the one-step ahead forecast errors computed from 

•t" Y fi - V l) 

where Y^d) is the one-step-ahead forecast from time origin t corresponding 
to the actual observation Y t+ ^. 
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Finally, tescs of significance with Q statistics reveal no systematic 
errors in the employment forecasts, suggesting that the observed errors might 
have been Induced by random disturbances. In contrast, test results do 
indicate that some systematic errors do exist in the Chinese Indicator 
forecasts. The associated Q statistic rejects the null hypothesis that the 
observed deviations between the forecasts and the actual observations are 
randomly Induced. 

Based on the above results, we must conclude that as of the date of 
the writing of this report there is still an absence of consistency in the 
forecast performance of the CEVMFC system. Further efforts are recommended 
and underway to improve the forecast accuracy of the system. 

A.5.3 Data Analysis 

Analysis of data with the CEUMFC system has produced some evidence in 
support of the hypotheses that 

(1) The frequency of demotion Increases during periods of political 
crises and displays patterns of collective experience shared by 
members of competing factions. 

(2) Changes in leadership resulting from factional differences lead to 
adjustments in policy priorities. 

Continuing the trends observed in the previous Rand study (Sung, 197S), 
the post-197A data reveal high demotion frequencies coinciding with known 
crisis periods. In particular, personnel demotion rate is highest in 1977 coin¬ 
ciding with the crisis attendant to the political demise of the "Gang of Four". 
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As cables 4a to 4e show, the total number of demotions in any one of the 
factional categories for that year is at least cvive as large as that for 
any ocher year. 

Moreover, significant differences also appear in the relative demotion 
frequencies between factions. High demotion frequencies appear to concentrate 
in particular factions. Among factions, defined according to functional 
affiliations, the highest demotion frequencies occur among military service 
corps (e.g., armor, artillery, infantry, public security, etc.) as compared to 
much lower frequencies among personnel in civil administrative functional 
areas. This reflects the continued effort to lessen military presence in 
the governmental and party hierarchy in the post-Lin Piao era. Among the field 
army factions, high purge frequencies are observed in the first, third, fourth, 
and fifth field armies. It is significant to note that the frequency of 
personnel loss is high in the 4th army with which Lin Piao was closely 
identified. In contrast, the 2nd field army faction continues to enjoy its 
post-Lin Piao era dominance observed in the earlier Rand study. It has the 
lowest demotion frequency of all the field armies. Likewise, a pattern 
of concentration also appears among the generational factions. Higher 
demotion rate seems to plague the older generations, such as generations one, 
two, and three, as contrasted to the relatively low frequencies in generations 
four and five. Thus, bevond personnel loss resulting from shifting power 
distribution, the foregoing results also show the slow but steady generational 
replacement effect at work. Moreover, civilian and military leaders have also 
fallen collectively as groups since 1974. With the exception of 1976, the power 
position of the military in relation to the civilian elites has undergone serious 
erosion since 1974. The pattern displayed is in line with the trend towards 
waning military Influence in party and government in the post Lin Piao era. 
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And consistent with earlier findings (Sung, 1975) no distinguishing patterns 
of behavior separate the relative demotion frequencies between commanders and 
commissars. 

Teats of hypothesis H., employing the models outlined In Appendix C have 
yielded some evidence of relationship between crisis occurrence and policy 
changes. A sample of the structural estimates for various models Is displayed 
In Tables 5a to 51. The stimates of crisis Impacts on policies are located In 
the verticle rectangular enclosures In the five series models and on the 
upper right-hand corner of the parameter matrices for the two series models, 
whereas estimates of the reverse effects are located In the horizontal 
rectangular enclosures in the five series models and on the lower-left hand 
corner of the parameter matrices of the two series models. Most of these 
estimates are significantly different from zero suggesting the existence of 
causal relationships between the crisis and the policy indicators. For 
example, the relative swings of policy preference between industrial and 
agricultural development seem to reflect the relative power positions of 
generation one and two and of civilian-military relations. Likewise, the 
patterns of external trade relations, measured in terms of the proportion of 
import and export to total GNP, also seem to respond to the sequence of crisis 
events. The crisis indicators, measured In terms of demotion frequencies by 
functions, by generation, by field army, by civil-military and by commander- 
commissar distinctions, all appear to exert strong Impacts on China’s level 
of foreign trade. Confirming some earlier findings (Liau, 1976) these results 
give further credence to the suspected Influence of domestic politics on China's 
relations with the external world. In short, our data reveal that political crises 
do have significant policy consequences In China. 
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Furchenaore, there Is also some limited evidence of a reciprocal 
effect of policy differences on inter-factional disputes. For instance, 
estimates in Tables 5a to 51 suggest that the ratios of industrial to 
agricultural outputs, of Import and export to GNP are causally antecedent to 
relative demotion frequencies In factions defined according to generational 
differences, functional differences, civilian-military distinction, and 
field army affiliations. These results suggest the possibility of policy 
differences fueling existing lnter-factlonal d putes. According to the 
power model this reciprocal effect would occur as factions highlight their 
policy differences to camouflage more basic differences over the distribution 
of power and political resources. Nevertheless, differences over policies, 
issue* and Ideology do serve to heighten the fundamental disputes over power 
dlstrubtlon. Yet, without more solid evidence, we can only tentatively 
entertain this power interpretation that political crises, which result in 
changes In leadership, carry significant policy Impacts and often contribute 
to policy changes, even though policy differences might have lnltally exacerbated 
the lnter-factlonal disputes that led to the outbreak of crisis. Uhlle a 
more precise Interpretation of our data must await further analysis, our 
findings have not been contradicted by the recent developments in China. 
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Table 4a 

DISTRIBUTION OF DEMOTIONS BY FUNCTION 







F u 

n 

c 

c 1 

0 

a 

n 





Year 

01 

02 

03 

04 

05 

06 

07 

08 

09 

10 

11 

12 

UNX TOT 

*1974 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

2 

1975 

2 

12 

1 

l 

0 

0 

0 

0 

2 

0 

0 

0 

7 

25 

1976 

l 

5 

0 

1 

l 

0 

2 

5 

3 

0 

1 

0 

2 

22 

1977 

15 

17 

0 

5 

1 

1 

0 

2 

2 

0 

0 

0 

38 

81 

1978 

5 

3 

l 

2 

0 

0 

0 

0 

2 

l 

0 

0 

15 

29 









Table 4 b 










DISTRIBUTION 

OF 

DDK) 

TIONS 

BY FIELD ARMY 


Year 

FAl 

FAl 

Fa3 


FA4 

FA5 

FA6 


UNK 


TOT 

*1974 


0 


0 

0 


l 


l 

0 


0 


2 

1975 


1 


1 

6 


5 


7 

4 


1 


25 

1976 


12 


1 

8 


10 


7 

2 


2 


41 

1977 


16 


9 

13 


13 


15 

12 


3 


81 

1978 


4 


6 

6 


5 


3 

2 


3 


29 


•Tha 1974 daca ara undarasciaaeas of the true demotion frequencies. 

A large number of demotion frequencies for 1974 were not considered In 
Che present sample. As the previous Rand study contained data up to 
September. 197} and since our present analysis began In January 1974, 
all cases In between were left out. Hence, until the missing data are 
accounted for our 1974 estimates must be treated with extreme caution. 







Table 4c 

DISTRIBUTION OF DEMOTIONS BY GENERATION 


Actual 

Year 

Cl 

G2 


C3 

G4 G5+ 

UNK 

TOT 

*1974 

1 

0 


1 

0 0 

0 

2 

1975 

6 

3 


8 

1 1 

6 

25 

1976 

3 

4 


7 

1 2 

24 

41 

1977 

11 

14 


15 

1 3 

39 

81 

1978 

3 

4 


6 

1 0 

15 

29 






Table 4d 





DISTRIBUTION OF 

DEMOTIONS BY COMMISSARS VS. 

COMMANDERS 

Actual 

Year 

CSR 


CDR 


BOTH UNK 

TOT 


*1974 

0 


1 


0 1 

2 


1975 

4 


11 


0 16 

25 


1976 

8 


6 


0 27 

41 


1977 

18 


18 


0 45 

81 


1978 

7 


8 


0 16 

29 




DISTRIBUTION 

OF 

Table 4e 

DEMOTIONS BY CIVILIAN VS. MILITARY 












































Table ’d 

Structural Model of Crisis and Policy Behavior 


Series 

Y i 

’t 2 

Yj t " lnduatrlal/agrlcultural 

2.28696 

.92357 

Y 2t - ailltary/GNP 

0.05087 

.00364 

Y 3t - Import s/GNP 

13.39696 

25.8210 

Y 4t - cxports/CNP 

14.20087 

21.046 

Yj t * army 3/army 4 

L.13609 

1.47215 


Final model - ARMAV(l.O) 


dependent 

Y lt-1 

Y 21— 1 

Independent 
l Y 31—1 

Y 4t-1 

Y 51- 

Y I 

' .69385 

.37355 

-.95428 1 

.0762 

-.032068 

Y 2 , 

. 3094 5 

.47646 

.99309 -1 

.2580 

.14412 

T 3 ! 

.57728 

-.06476 

.28414 

.25812 

.077751 

y 4 j 

. 384 08 

. 35024 

-.14481 

.88775 

.026938 

Y 5 I 

|-.91837 

.34083 

.37253 

.10855 

-.16016 
















Table Sf 

Structural Model of Crisis and Policy Behavior 


58995 


12550 


.84583 


r 

Series 

*1 

1 

.ill. 

Yj t - Industrial/agricultural 

2.2869 

.92357 

Y 2t - milltary/GNP 

.05087 

.00364 

Yjt - laporta/CNP 

13.396 

25.2821 

Y^t ■ exports/CNP 

14.2008 

21.046 

I Yjt - civillan/milltary j 

1.55652 

4.83198 
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Table 5g 

Structural Model of Crista and Policy Behavior 




2 

Series 

Y i 

°i 

Y^ c “ Industrial/agriculture 

2.28696 

.92357 

Yj t • army 3/army 4 

1.13609 

1.47215 


Final 


Y 

Y 


It 

2t 


Y 

Y 


It 

2t 


Y 

Y 


It 

2t 


Model - ARMAV(2, 

Y lt-1 
.10862 
. 54098 

Y lt-2 

:~7l8493 

2J3590 

a lt-l 

[-.*87145 

I 

. 2.836 


1 ) 


-.43470 
-. 15756_ 

Y 2t-2 

-. 14847” 
-. 11882_ 

a 2t-l 
-.53762 
-. 58577 
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Table 51 

Structural Model of Crisis and Policy Behavior 


Series 

? 1 

2 

° 1 

Yj c " induatry/agriculture 

Y 2t " civillan/ailltary 

2.28696 

1.55652 

.92357 

4.83198 


Final Model - ARMAV(2,0) 
Y lt-1 

Y lt .84853 

Y,, .73531 


lt-1 

.065049 


.01413? 

. 6684J__ 

Y 2t- 

.080411 

-.46098 





































A.1.0 Autoregrcaalve-Hoving Average Vector Time Series Model 

The "autoregressive-moving average vector" (ARMAV) model In this software 
describes a system of feedback relations within the framework of an autoregressive- 
moving average process. All variables are assumed to be interdependent over time, 
so that the system encompasses a range of structural relationships. Nonetheless, 
a parsimonious description of the system is possible through vector representation. 

Consider a system in which the state variables at time t are dependent on 
variables at preceding times, t-m through t-1, and on their corresponding noise 
series from times t-n to t-1. Then a representative autoregressive-moving average 


(A.1-1) 
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s^ “ seasonal operator of Indicator 1 

■ number of seasonal differences for Indicator 1 
dp - number of regular differences for Indicator 1 
The subscripts In the above equations denote the row and column position 
of elements In the system of equations and the order of the autoregressive or 
moving average process respectively. In each of the above matrices, ♦ or G, the 

diagonal elements. Qr (for 1-1.p. j-1.m and k-1... .n) contain 

the autoregressive and moving average parameters for each series, l.e., 


th. <Up.nd.oc. of Y u on T lt _,. of T ( on T u ... « Uo of T t on 

and ? lln ° f Y t 00 *1t-n' Th * off-dlagonal elements. In contrast, determine the 
Inter-relationships between variables. For Instance t m defines the dependence 

on ^2t-l* *12m t *' e dependence of on ^2t-m *211 the dependence °f Y^ on 

T lt-1* * 21m ° f T 2t ° n Y lt-n - A similar Interpretation applies to the off-diagonal 
elements In the 0 matrices. Elements In the (I - B) * matrices are difference 


operators (l.e. (l-B)Y t » Y £ - T^^). Use of these onerators In time series models 
Is one wav of ellmlnatins nonstationary In time series. Svstems so differenced will 
have characteristic roots that lie within the unit circle. The condition for 
statlonarlty In the vector model will be that the eigenvalues of the : and 0 matrices 
lie within the unit circle. This condition will be periodically checked in our 
modeling procedure. 
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Each element in ^ is assumed to be normally distributed with mean zero 

2 

and variance Moreover, ^ also satisfies the conditions that 

(1) is Independent of g £ _ n 

(A.1-2) 

(2) ^ U Independent of . 

In other words, no reverse dependence of ^ on past ^'s or on past ^'s beyond 
the autoregressive order is permitted. Furthermore, since ^'s coefficient matrix 
is the identity, no contenponneous relationships are assumed to exist smong the 
elements of ^ (e.g., ••• •Xp t )• Finally, no restriction is Imposed on the 


’it 

Moreover, by imposing certain restrictions on elements of the vector one might 
also entertain various alternative special cases. For example, the existing system 
is reducible to a single equation autoregressive-moving average process: 

Y t + *m Y t-m * *t °l*t-l“ ~ Q n*t-n (A.1-3) 

Although other models are reducible from (3.1-3), further analysis of their 
likelihood function is necessary before we may implement them In our forecasting 


software. 
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A.2.0 Modeling Strategy: Fitting Autoregressive-Moving Average Vector Model 

Our modeling strategy encompasses a three-step procedure. The steps are 
identification, estimation, and diagnostic checking. In identification, we 
select that model which best fits the observed time series. In estimation, we 
derive estimates of the model parameters. In diagnostic checking, we check the 
adequacy of the estimated model. Inadequate models are re-identifled and re- 
estimated until adequacy is assured. 

Moreover, ldentiflcation, estimation, and diagnostic checking are performed 
within a hypothesis testing framework. By setting up pairs of successive orders 
of models under null and alternative hypotheses, we can test alternative choices 
of parameter restrictions. Through successive comparison of higher order models 
we will be able to Incrementally approximate the observed time series and determine 
the best combination of autoregressive-moving average forms. 

A . 2,1 Model Identification 

Employing a hypothesis testing approach, vector model identification calls 
for the fitting anc comparison of pairs of successively higher orders of models 
until all of the systematic variations in the time series are accounted for. 
Accordingly, successive orders of ARMAV (m,n) models are fitted until a point 
is reached beyond which no statistically significant improvement in model 
approximation can be expected. The last model to register a significant improve¬ 
ment is retained as the optimal or best fitting model. This procedure is outlined 
in the following steps. 

(1) The first question to be posed in model fitting is: what is the optimal 
order of the ARMAV (m,n) model? In other words, what are the values of m and n? 
This question is best answered by examining successive values of m and n until 
a model is reached that no longer displays significant Improvement in goodness 
of fit (this criterion will be defined later). This approach presupposes a 
systematic and natural progression in the sizes of m and n examined. Indeed, 
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lc has been shown chat any stationary process can be adequately approximated 
by Increasing orders of an ARMAV (m,m-l) model for large enough m (Pandit, 1973). 
Accordingly we Introduce a modeling strategy based on the following sequence 
of steps: 

a. Beginning with ARMAV(1,0) we Increment our model by (m,m-l) order, 
where 01-1,2,.... Along this progression, we would successively examine 
higher order models, such as ARMAV(l.O), ARMAV(2,1), ARMAV(3,2), ARMAV 
(4,3),..., ARMAV(m,m-1). The Increments are by order of 1 for both 
Che autoregressive and the moving average components. Ue terminate 
this progression when no further Improvement in goodness of fit can 

be expected from proceeding to a higher order model. 

b. Since the models are Incremented by orders of (m,m-l), there are 

alternative Intermediary models not accounted for by the above 
procedure. These are models resulting from Incrementing the m or n 
value separately by order of 1. These alternative possibilities to 
each ARMAV (m,sr-l) model are shown In the boxes enclosed In dotted 
lines as shown In Figure 3. Within each box we compare the 

alternative Intermediary models to the mainline ARMAV (m,m-l) models 
.to determine the optimal model according to some statistical criteria 
of goodness of fit. However, we do not examine all the boxes, but 
only the box corresponding to the best ARMAV (m,m-l) model. Thus 

we begin by successively incrementing the autoregressive and the 

moving average components In ARMAV (m,m-l) model by orders of 1 

until we arrive at a mainline model beyond which no further Improvement 

can be expected. We then test the alternative Intermediary models In 

the corresponding box against the final ARMAV (m,m-l) model, and determine 

the optimal choice. For instance, we might first proceed along the 
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mainline models (e.g., ARMAV(1,0), ARMAV(2,1), ARMAV(3,2), AKMAV(4,3)) 
and find that the no significant improvement In goodness of fit results 
from going beyond AIMAV(3,2) to ARKAV(4,3). We would thus drop back 
to ARMAV(3,2) and test It against its corresponding intermediary models 
enclosed in the same box. By focusing on the Intermediary models in 
the same box, namely, ARMAV(3,3), ARMAV(4,2), ARMAV(2,2), and ARMAV(3,1), 
we have successfully bypassed suboptlmal lower order intermediary models. 
The saving in terms of computational effort is obvious, 
c. As a final test of model adequacy, we check for superfluous moving 

average parameters in the final model selected in the preceedlng step. 

In this test, we successively decrease the order of the moving average 
parameters until no superfluous parameters remain, as a substantial 
change in model goodness of fit, due to the dropping of parameters, would 
indicate. Accordingly, the final check Involves testing in declining 
order a sequence of models: ARMAV(m,n-2), ARMAV(m,n-3),..., ARMAV(m,0). 
(Note that this final step is not necessary if the optimal model selected 
at step b is ARMAV(m.m-l), since its corresponding lower order moving 
average models also constitute the intermediary cases already examined). 
Methodologically, this approach adopts a regression perspective. It is 
predicated on the assumption that the dependency in the time series can principally 
be modeled by autoregressive processes. The moving average processes merely account 
for the remaining stochastic dynamics in the series. Thus, by this reasoning it 
is entirely possible to couple a lower order vector moving average with an ARV(m) 


component 




ARMAV(m.n) Modeling Strate 


ARMAV(m,m-2)-—>> 


ARMAV(4,4) ARMAV(5,3) —'“H —> 

' ^ S’ ' 

I ___ ARMAV(4,3) I 

|_ P—ARMAVO, 3)^1. ARMAvdT:!”-^^ —> 

|_\\RMAVO,_2_)_ | J 

l _l_ARMAV(2 1 2)^_ / j^_ _ AKMAVO.T? i -> ft 

' I jn ? 

^ARMAV(2,1>\,_ j | 

j^ARMAVU.l)^- _ A _ _ ARMAV(2,0)j-> tj 

r i 

ARMAV(l.O) 



check for superfluous 
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(2) Selection or model is based on a goodness of fit criterion (smallest 
sums-of-aquares of the residual) by which we Judge the adequacy of a model in 
capturing the systematic trends in the time series. The goodness of fit criterion 
is applied through a hypothesis testing approach. Accordingly, we examine a pair 
of models at a time in increasing order to compare the observed difference in 
their goodness of fit and to determine whether the differences are results of 
random errors or of some systematic effects. The higher order model is known as 
the full model and the lower order model as the restricted model. Parameters 
which appear in the full model but not in the restricted model are assumed to 
have been restricted. Applying the hypothesis testing approach, we would set 
up a null hypothesis to correspond to the restricted model and an alternative 
hypothesis to correspond to the full model. The null hypothesis defines the 
restricted parameters as zeros and the alternative hypothesis sets them to some 
non-zero values. Thus in comparing a restricted ARMAV(2,1) model. 


to a full ARMAV(3,2) model. 


(A. 2-1) 

(A.2-2) 



(3) A statistical test of significance serves to distinguish between the 
full and the restricted model. The choice of model is based on two alternative 
tests. They are the F test and a test involving comparison of forecast mean- 
square-errors. We describe these tests briefly below. 

(a) The F test determines whether the observed differences in the models' 
residual si»-of-squares, which measures the models' goodness of fit, 
indicate true Improvement in model adeouacv or are attributable to 
random effects. The F statistic used is (Rao, 1973, pp. 555) 
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F Testa for Alternative p and a 


Values of p.s 

F ratio test 

degrees of freedom 

s-1 for any p 

(1 - A) /p 

p and N-r+s-p 


A/ (N-r+s-p) 

(A.2- 

s-2 for any p 

Cl - ''fi/p 

2p and 2(N-r+s-p-l) 


/\7 (N-r+s-p-l) 



where 

r - number of parameters in the full AR.MAV(ra,n) model - (m * n)p 2 
s - number of parameter matrices restricted to zero in the restricted 
model 

p - number of variables or series in the vector model 
N - number of observations It, each time series. 



Th e |\J and| A A |are determinants of the quadratic residual sum-of-squares-and- 
product matrices (i.e., • a^a^ where a is an nxl vector of noise series 

containing elements for the restricted model (under the null 

hypothesis) and the full model (under the alternative hypothesis) respectively: 
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b. Forecast mean-square-error comparison permits us to determine Che 

relative forecast performance among any number of models. According 
to this test, the model with the least forecast mean-square-error 
has the best forecast performance. We use for this test a statistic 
known as the "Thiel coefficient" (Thiel, 1965): 

MSE - l 2 + + S y - 2 S y S y r(Y t .y t ) 

where Y f is the time series observations and y £ 
corresponding forecasts. 




are the square of the averaged forecast errors, the standard deviation of , 
the standard deviation of y^t and the correlation coefficient between Y^ and y^ 


(A.2-5) 

- Y <t) their 


respectively. 
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By this test criterion, we would select as the optimal model the one giving 
the least forecast mean-square-error measured In terms of the Thiel coefficient. 

(4) The foregoing tests are used alternately In different phases of our 
modeling strategy. Within the hypothesis testing framework the F statistic Is 
used initially In selecting the best mainline ARMAV(m,m-l) models and then is 
used in identifying the "significant" Intermediary models localized about the 
best ARMAV(m.m-L) model within a given box in Figure 3. Having Identified the 
significant models we next compare the estimated Thiel coefficients for each 
model to determine the model which gives the best forecast results. Finally, 
we test for superfluous moving average parameters in the optimal model by 
comparing its forecast mean-square-error with those of models having successively 
lower order moving average components. This procedure continues until no 
superfluous parameters are left. The model that we ultimately retain is 
therefore both significant and has the Least amount of forecast errors. 

A.2.2 Model Estimation 

Model estimation entails an automated three-step procedure. In step 
one, model fitting is accomplished by linear least-squares. The CEWMFC 
system automatically fits data to successively higher orders of autoregressive 
models. In step two, a transformation, through an "inverse function", of the 
autoregressive parameters generates initial estimates of nonlinear autoregressive- 
moving average vector model parameter. These estimates are needed to start off 
the iterations in nonlinear least-squares estimation. In step three, final 
parameter estimates are derived for the autorerressive-moving average vector 
model. Using the estimates from the inverse function as initial values, a 
nonlinear least-squares routine, based on a FIetcher-Power minimization 
algorithm, fine tunes the initial values until they are near maximum likelihood. 
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A.2.2. I Estimating ARV(m) Linear System (ARVEST) 

Estimation of a pure autoregressive vector model may be accomplished by a 
linear least-squares method. Consider, for instance, the model 


w h«* r « X t , and a t are lxp vectors; and $ is a pxp matrix of coefficients. To 
estimate the { matrix we minimize the sum of squares and cross-product matrix for 


A - E (Y t - ♦X t )' <Y t - ? X t ) 


Differentiating A with respect to {yields. 


Setting the derivative equal to zero and substituting i for 


get the normal 


from which the least-squares estimate obtained is 


In estimating the vector model 

*0 ?t - *l?t-l + ••• + fJt-m (A.2-11) 

we set X t as a partitioned matrix: 

- [it; Vi -t-mi 

and regress T on X as shown above. Setting we get 

w *c -1 ** I>U -m -m*0 

the linear ARV(m) estimates. 
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A.2.2.2 Est imating Initial Values for ARMAV(M.g) Itode1 (INVEST) 
the initial values for the ARHAV(n.n) nodal may ba derived by making use of a 
simple relationship between a mixed autoregreealve-moving average process and 
a pure autoregressive process. It is known (Box and Jenkins, 1970) that a mixed 
ARMAV(m,n) process, expressed in terms of backshift operators, B ■ (B^.B p ). 

(X - fiB - ... - • j B*)Y t • (I • B • ••• “?n? n, -t CA.2.12) 


has an equivalent "Inverse" function. 



where the *'s where the linear combinations of the ♦'s and 9's. 

The above relation Is established by substituting the pure model for a ( 
mixed model, which gives us an operator Identity (after cancelling out 
both sides of the equation): 


(A.2.13) 

t In the 

T t fro * 


(l - ?l i - ... - ?- B S ) - (I - ?jl - ... - •/*> (I - Tl B - !2 B 2 - 


...) 


Equating coefficients of equal powers of B (see LI, 
we get the following system of equalities 
!i 

U • ?2 - *iV !2 


1979; Box and Jenkins, 1970) 


(A.2.14) 


-i " ? J ’ " -2-J-2 - ♦ *J 

In particular, for ) greater than either of the larger orders of a and n, we 

have 

<! - !i* - M*.!.**> m 

It should be apparent from the above expressions that the Initial guess values 
for the f's and the «’s can be derived once the *'s are known. Since the Inverse 
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function Is an Infinite order autoregressive process, it can be 
approximated by a vector autoregressive process: 


- (I - i.B- ... - r B )Y 


and defining k ■ m+n, we obtain our estimates for the Inverse function. This 
equality holds so long as the moving average parameters are invertible, that is. 
If the moving average component converges. This requires that the eigenvalues 
of the •» matrices lie between -1 and 1. When this condition is not satisfied 
by the data at hand we may impose convergence on the moving average vector 
component by inverting the eigenvalues that are larger than 1,1.e. i/c^. to 
assure compliance with the invertibllity condition. Since the above vector 
model is linear we may estimate the r'a by a linear least squares method 
£ - <X’X/ 1 X*T 

C C C ' (A.2.17) 

where X - [T ! Y • . . . Iy 


derived our initial estimates for the j’s. 

Consider, for example, an ARMAV(2,1) model: 

-t * -l^t-1 * * 2 * 1-2 * *t " ?l-t-l (A.2.18) 

In order to derive the initial estimates for and 9^, we must first approxln 

the corresponding inverse function. Since there are three matrices of parameters 
in the above equation, the inverse function would be 

?t • -i-t-i * !2 T t-2 ♦ ?J?t-3 
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for k - 2+1 - 3. Using linear lease squares method we derive the linear approx¬ 
imation for B^, S,, and 8^. Then, equating and using (A.2-15) for 

J > 2. we derive our Initial estimate of as follows, 

(I - e^) Tj ■ 0 • for ) * 2 

which Is 

(I - ^B) * 3 -0 



Solving the above equation for 4^ we get 



Next, we check for Invert lb11lty. Invertlblllty requires that the eigenvalues 
of lie within the unit circle, -1 and 1. If any eigenvalue exceeds the 
upper and lower limits we would Invert It. (e.g., c* • 1/*^). We compute the 
eigenvalues using the RTFIND routine. We substitute the c* for the Invertible 
eigenvalues and reintroduce the results Into the following expression 

- **! (A.2-20) 

where T Is the pxp matrix of eigenvectors corresponding to the original untransformed 
'lgenvalues, and • * is a pxp diagonal matrix containing inverted eigenvalues. ©* Is 
the new coefficient matrix to be estimated from T and c*. Transposing matrices 
In the above expression we may obtain the new matrix: 

!f • T**T” 1 

Since c contains only invertible elements would satisfy the lnvertiblllty 
condition. This step will be applied to the matrix even when it contains no 


invertible eigenvalues. 
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Knowing *2» a nd we can then calculate and 

Using (3.2-12) we 8«t 



Accordingly, for an ARMAV(2,1) model containing two enuatlons, let the least squares 
estimates from an ARV(3) model be 



Finding the eigenvalues of 
V 2 - . 898X «• .258 - 0 



Note that 3j i« ass-imed to be equal to tero in the last equation. 
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Node* Chat Che eigenvalues of 9j, are {-.572, .466}, with magnitudes less 



.3 Estimating ARMAV(m.n) Sonllnear System (SMLEST) 


Estimation of a nonlinear AKMAV(m.n) model Is achieved by minimizing the 


determinant of the sum-of-squares-and-product matrix of a where 


In which Y^'s and a^'s are pxl vectors; the $'s and the 9’s are pxp matrices 
of autoregressive and moving average parameters. The matrix whose determinant 
Is to be minimized Is 



St • 1/n A 


(A.2-24) 
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The method for minimizing A with respect to the ♦ and <3 parameters used 
in this software system is known as the "Fletcher-Powell" method. It is a 
modified gradient search routine, which, in contrast to the popular steepest 
descent method that makes use of first derivative information, utilizes second 
partial derivatives to determine the most direct path by which to move some 
initial estimates (from the INVEST routine described in the preceding section) 
towards their maximum likelihood values. 
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(6) We then calculate the second derivatives from differences of the gradients 
for evo successive iterations: 


(7) Using the iC and 5 vectors computed above we next construct the A and 1 


(8) Finally, we compute the Hessian matrix (note that H Is an Identity matrix 
before the first Iteration) according to 


(A.2-37) 

(9) The above procedure Is repeated for each succeeding Iteration, during which 

the H matrix Is successively approximated. Iteration continues until changes 
In the difference vector and the direction vector between two successive 
Interatlona are less than or equal to some predetermined level, c, for In¬ 


entire procedure replacing for the previous H g . 

Since minimizing the differences between the AG g and Is equivalent to 

objective function over a suitable choice of t and 0 parameter 
values, the above procedure will lead to maximum likelihood estimations of 


the autoregresslve-movlng average vector model. 
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.K.2.2.U Minimization Criterion 

The determinant of the sum-of-squares-and-products matrix of a t> A^a^a^, 

where a - (a, , a. . a 3, the vector of residuals corresponding to 

s It 2t pt 

equations 1 through p. Is 

|-*lt U lt*2t •" E *ltV I 


pt It pt 2t * pt 


which by denoting each entry as d,. for p and J“l,...,p Is 


d ?l d 22 d 2o 


LI ’ 22 pp 


(d_ .)«._> 


In a two by two matrix (for ARV(l) model with p-2) the determinant of the squares 
and product matrix to be minimized is 


It 2t 'It 2t 
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which oust be positive if Ag Is positive definite, as Ia^ t a“ t is alwa y s greater 
than E(a lt a 2t ) 2 . (Coldberger, 1964) Minimization is achieved by finding that 
set of estimates giving the smallest D^. 

A.2.2.S Matrix Inversion and Characteristic Roots 

Several subroutines are designed to compute the determinant and the inverse 
of a matrix and to compute the characteristic roots of the ARMAV(m,n) model in 
terms of the eigenvalues(and eigenvectors) of the parameter matrices. The matrix 
inversion routine is call"MATNVl" and the root finder is call "RTFIND". They are 
used in linear-least-squares estimation, initial value estimation (INVEST), and 
nonlinear estimation among others. The determinant of a matrix can be numerically 
estimated by an expansion: 


a, C, 
ip ip 


(A.2-40) 


' ;.,vu 


where a^ is the element of the matrix in the ith row Jth column; is its 
corresponding cofactor which is defined by the relation 




and | is Che determinant of the minor of Che cofactor. The minor is the 
submatrix that results if we eliminate the ith row and the Jth column (inter¬ 
secting at a^) of the matrix. 

The inverse of the matrix accordingly is simply 

A" 1 - C’/I^l (A.2-41) 

where C' is the transpose of the cofactor C of matrix A. For such an inverse 
to exist would require that A be a square matrix. 

The eigenvalues and eigenvectors to the A matrix are solutions to the 


following system of homogeneous equations 
(A - XI)T • 0 


(A.2-42) 
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where X Is a lxp vector of eigenvalues and T is a pxp matrix of Lxp eigenvectors. 
Nontrivial solutions for the above equations exist if the determinant 
of the matrix of coefficients (A - XI) is equal to zero: 


|A - Al| - 0 


(A.2-43) 


In general, the eigenvalues can be real and distinct, real but equal, or 
complex. Real roots define dynamics displaying the sum of exponential patterns 
of behavior; complex roots define cyclical patterns of behavior. For a two 

equation ARV(l) system, for example, the roots or eigenvalues are 

/ 2 ’ 1 (A.2-44) 

. . *111 ’*221 /("ill " *221) *221 *121 

V*2 2 2 


These roots are real if the "discriminant" under the square root is positive or 
equal to zero: 

2 

“ ( * 1 U ‘ * 221 ) *211 *121 - - 0 

and they are complex if the discriminant is negative: 

C< *111 - * 221 >'-‘ *221 * 121 ^ ’ 

When roots are complex they form complex conjugate pairs such as 

i i ^ ♦! (A.2-45) 

* 1 * *2 • d - 

- d (cos 8 ♦ 1 sin 8) 

where 1 - J-T, d • Xjj ♦ X^, 8 - tan’ 1 (——). The subscripts I and R stand 
for the real and the imaginary parts of roots into %4Uch complex roots may be 
split, for Instance, the complex conjugates X^ and X^ may be written as 

X l* X 1 " X R - 1 X 1 (A.2-46) 

For further discussion see Li (1979, ch. 10). 
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A.2.3 Diagnostic Checking 

Various diagnostic statistics are Included in this system for testing the 
goodness >f fit of the estimated ARMAV(m.n) model. First, Inadequacy of the 
autoregressive component of the model Is determined by examing the residual 
autocorrelations and partial autocorrelations for each equation. The autocorrelations 
are computed according to 


I (a - a ) (a 
„ , t-1 

r (It) . - 


The partial autocorrelacions are recursively computed from 


A two standard error confidence limit is defined for the abo< 


alues exceeding this limit suggest inadequacy of the model. 

Moreover, since it is known that the residual autocorrelations are inefficient 


timates and tend to understate 


at the lower order lags, an overall test 


of significance is also introduced (Box and Jenkins, 1970). The statistic is 


which is distributed as chi-square with k-m-n degrees of freedom, where k is the 
number of lagged correlations, m the order of the autoregressive component and n 
the order of the moving average component and N the number of observations. 

Second, any needed revision in the interdependent relationships between 
variables or series will be revealed by examining the cross-correlations between the 
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residuals and the input series. The sample cross-correlations are computed 
according to 


l (a 

k-0 




{“ 


(a - a r E (Y 
k-0 C C k-0 C 


A two standard error confidence limit for the above is 

r ay (k) + 2 \j~l/ (n-k) ’ (A.2-52) 

Values exceeding this limit will again suggest model inadequacy. 

An overall test across all cross-correlations is also provided in the 
following statistic 

s “ N ^ r av < k > (A.2-53) 

k-1 ' 

which is distributed as chi-squares with k+l-(mp‘'-l) degrees of freedom, where 
k is the number of lagged cross-correlations and (mp 4- -!) is the total number of 
parameters in the autoregressive component of each equation in the system. Model 
inadequacy is suggested when the S statistic is larger than a certain critical 
chi-square value determined by the routine to be described in the ensuing section. 
A.2.3.1 Chi-square Statistic 

In order to determine the significance of the Q and S statistics we use a 
chi-square table to find the critical value up to 30 degrees of freedom, at 5X 
significance level, and for greater degrees of freedom we use a standard normal 
approximation of chi-squares: 

S ( Z 3 ♦ V 2d-1) 2 (A.2-54) 

where d is the degrees of freedom from the Q and S statistics and a - 1.645 at 

5X significance level. The search routine for determining the critical chi-square 


values la programmed in the CHIS<) routine. 



A.3.0 Forecasting 


The Forecasting Routines eaploy the following algarlthms. 

A.3.I PSI Weights (PSICOF) 

(I - ?l B - } 2 B 2 - ... - fjfti * + *B + , 2 B 2 + ...)- (I - 0 r B - ... - e n B°) (A.3-1) 
Equating llke-terms In the above equality gives the following system of equations: 

*0 " ? 

" Mn " 


:i -i-o -i 


*: *o " h 


-2 * 1 J 1 *2 y -0 


*3 " *1*2 + *2 *1 * *3*0 ‘ : 


*»-l " tl*-2 * *2**0 + + *m-l*0 

+ ♦,*, ♦ ... ♦ ♦ * n - e 


-n-l 


-a -1-m-l -2-m-2 *m-0 -n , for m-n 

*m " !l*m-l + *2*.-2 + ‘ • * + *m*0 • for a>n 

Setting the Initial PSICOF matrix equal to Identity, we aay recursively compute 
higher order matrices, knowing In advance the estimated } matrices. For example, 
in an ARMAV(2,1) model where the estimates are tpfj* and Qj* th « PSICOF expansion Is 
*0 " - 1 


-1 

*2 " * 1*1 


1*0 


- 2-0 


|lVl ♦ *2*.-2 


So, for 
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These matrices are used In several parts of this software. Specifically, they are 
necessary In computing (1) the forecast confidence Intervals, and (2) the adaptive 
forecasting routine. 

A.3.2 Forecasting Algorithm (FORCST ) 

The forecasting algorithm Incorporates two major features. First, It permits 
forecasting with leading Indicators using multiple series. The forecasts are 
generated either from the full vector model or from a restricted single 
equation model. Second, It accommodates regular and seasonal difference operations 
to restore statlonarlty to the original scries. The difference operators are 
explicitly Incorporated In the forecasting functions. 







Mathematically, our problem is one of expending a forecasting system of the 




The general system forecasting equation Is: 

p • total number of series used In forecast 1 < p < 5 
*ljk " parameter of the k th matrix. In row 1, col J 

®ljk ” ■ ovin 8 average parameter of the k th matrix, in row 1, col J 
B • backspace operator; BY t • 


teasonal operator of Indicator 1 












Using the matrix column portion 

d., 


' c .»u I Ui 1 -• 1 Si 3 • c * lll ”» \t ■ 


: 31 I - I O 


h Y lt 

' ♦l^lt-l 

" *12 Y lt-2 * " 

' * + *ln Y lt-n 

X 2 Y 2t 

* *21*21-1 

' *22 Y 2t-2 + ’ 

• ^2n Y 2t-n 

I Y 

- » ,Y 

- » „Y „ - .. 

.. - a Y 

P Pt 

pi pt-1 

p2 pc-2 

pn pt-n 


- 3 - 4 ) 


Note that this expression of the system breaks the system into p equations, where 

each equation i depends only upon the i column of each coefficient matrix 

and the ith indicator. The differencing within each indicator can then be expressed 


by expanding the vectors $ in the manner Indicated bv the expression 

# i d il *12 

(1-B l ) 11 (1-B) . 


T1 
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1 * *111 " *112 
*211 + *212 


*121 

Y it-3 : 

- 

* l **222 

b-J 



+ *112 
+ *212 


*121 
1 **221 




•*121 * *12z" 


- 

u 

'*122 

"^2t-6 

”*221 * *222 

Y 2c-5 


0 

"*222 

Y 2t-6 


A.3.3 Forecast Confidence Intervals (FRCONF) 

A two standard error confidence Interval is cstiaated for each forecast. The 
confidence Interval la: 

T t (D i 2 /D^TTT 1 (A.3-5) 

where D^ t (t) Is the 1th diagonal element of the variance-covariance matrix of 


forecast errors for t-step-ahead forecasts: 


(A. 3-6) 


The w^’s are the PSICOF (see preceedlng section) and j* can be approximated by 
the sample residual sum-of-squares-and-products matrix. 


U 2t*lt U 2i 


It pt 
*2t a pt 


pt It “ pt 2t 


A.3.4 Forecast Updating (FRUPDT ) 

In many operational settings in which forecasts are made on a continuous 
basis. It Is expedient to update forecasts from time to time with the latest 
piece of information without recomputing the previous forecasts. This can be 
readily done for an autoregressive-moving average vector model, requiring only 
the previous forecasts and the new updating information. The method Is known 
as adaptive forecasting. * 
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Adaptive forecasting seeks to improve the accuracy of present forecasts 
by means of the current one-step-ahead forecast error computed from the latest 
updating Information. Assuming that we have a series of t-step-ahead forecasts, 
generated from time t, and a new piece of data becomes available at time t+1; 
our one-step-ahead forecast error would be Y t+ ^ - Y^d), the difference between 
the new data from time t+1 and the one-step-ahead forecast for time t+1 from time 
t. By properly weighting this measure of forecast Inaccuracy and adding the 
result to the one-step-ahead forecasts we may obtain significant improvement over the 
previous forecasts (Box and Jenkins, 1970). Thus, our adaptive forecasting 
algorithm for our vector model is 

' !t (t ♦ l > * ' !t u):l 

where Y c+ ^(l) is a pxp matrix of updated forecasts for time t+1; Y t (i +1) is a 

pxt + 1 matrix of previous forecasts generated at time t; Y^ + ^ - Y t (l) is a pxl 
vector of one-step-ahead forecast errors; is the PSXGOF matrix for the 1th 
step ahead forecast, containing nonzero diagonal elements and zero off-diagonal 
elements. 

The practicality of this method stems from the fact that in order to make 
use of the method we need only to keep track of our previous forecasts and of our 
one-atep-ahead forecast errors, albeit we might periodically want to reexamine our 
forecasting models with the original data. In short, this method permits us to 
learn from past errors in improving our present forecasts. 

A.3.5 ZStatistic and Standard Normal Probability Function (PRB2 ) 

A standard normal "Z" statistic is computed for each forecast. It provides 
a standardized measure of the devlstlon of a random variable from its mean. Taking 
the mean of an indicator or a series to imply its norm, the Z statistic may be 
construed as an indication of the degree of abnormality, or deviation from norm. 


in our forecasts. Assuming statlonarlty for our series, we may expect the 
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and che variance of our aeries Co remain constant over time. Thus, even though 
the mean and the variance of the future values that we are attempting to forecast 
are unknown, we may use the existing mean and variance as approximations. 
Accordingly, Che Z value corresponding to a forecast may be computed from 


where and are the sample mean and sample standard deviation of the existing 
series, and Y £ (t) Is the Ith-step-ahead forecast generated from time origin t. 
Statistically, Z is known to be normally distributed with mean zero and 

variance one. Its probability density function Is 

1 -z 2 /2 

f(Z) " ^ e * (A. 3-10) 

Integrating this function from 0 to Z would give us the probability for the region 
between the mean and the given Z^ value. The probability that a deviation from 
the mean will exceed Z^ Is given by subtracting the integral for the region between 
0 and Z from .3; chat is 


v - ,ir 


In the CEWM forecaster, the computation of the Z value Is performed In the 
main program and of the probability estimate is computed In the PRBZ routine. 
Subroutine PRBZ performs the following computation. 

The ctaulatlve density function 

F( X) - 7 J-* / X • dv (A.3-12) 


is. In general, difficult to evaluate for an arbitrary value x- The approach 
used In this system Is to estimate the value of F(x) by creating an Interpolating 
polynomial. By picking points (x^t^). i • 0, n and Y^ - F(x^)t and fitting 
a polynomial P (x) through these points, such that P (X,) * F(x.) a good estimate 








Using the data points In Che cable be'.ov (F (x) values are from Che 21st edition 
of the CRC Standard Mathematical Tables): 


1 *1 Y i - f <*i> 

0 

0 

.5 

1 

1 

.8413 

2 

2 

.9772 

3 

3 

.9987 

4 

4 

1.000 

5 

.5 

.6915 

6 

2.5 

.9938 


F^(x) UAS found to be 

P 6 (x) " P 5 (x) + *00050608 ( x ) (x-1) (X‘2)( X -3)( X -4)( X -5) 

- Pj(x) ♦ .00050608 X (X 5 - 10.5x 4 + 40 X 3 - 67.5 X 2 ♦ 49x-12) 

- .5 ♦ .390295221x+.0412013U x 2 - .137276229 X 3 
♦ .055993667 X * - .00951905l x 5 + .00060508 x 6 

Empirical coaparIsons have shown that this polynomial has a maximum error of 
approximately . 001 , which was considered to be sufficiently accurate for our 


purposes. 











A.3.6 Combined Forecasts (COMBFR ) 


When several alternative forecasting methods are available it is sometimes 
useful to combine their results in order to complement their strength and 
weaknesses. The forecasts to be combined must come from the same series though 
they are generated by different forecasting methods. Thus, the method is 
intended to compensate for the differences between methods, not between forecast 
Indicators. 

The method we use, as originally developed by Granger and Newbold (1969, 
1974), permits us to combine multiple frequency forecasts. The combined 
forecast is a linear combination of the original forecasts: 

C .t " *tV° (A. 3-14) 

where, K - (k, , k, . k ), is the combine forecast weights, and Y (£) 

it zt nt t 

is a matrix of alternative forecast series. 


Granger and Newbold (1974) suggested chat the forecast weights oa 
computed according to the following relation 

T '‘ 2 ,-1 , 

’it > 1 (or l-l, 


T-l 

. • ( l 

t“T-v 


2 i-l 


) 1 / ( r ( i 

1-1 t-T-v 


(A.3- 


where e - Y -Y (l), which is the forecast errors for £-step-ahead forecasts, 
it it-* it 

It is suggested that the forecast error for the last twelve periods be used to 

compute k it , which means in the above equation we set v-12. Our decision rule 
accordingly specifies that if the number of our series from which we generate 
the forecasts has more than 90 observations we set v-12; otherwise we set v-.12n, 
where n is the number of total observations in the time series. 














The crisis Indicators are developed from data on 17 biographical 
attributes. They Include: 

- name 

- date of death 

- date of birth 

- place of birth (native place) 

- generation 

- field army affiliation by generation 

- coonander or commissar by generation 

- date entered party 

- military or civilian 

- combat experience 

- awards 

- civilian education 

- military education 

- military-region affiliation by generation 

- functional affiliation by generation 

- position level or rank by year 

- province or military-region affiliation, by year. 

These attributes are selected according to criteria pertinent to data 
availability, institutional loyalty, and Interest group behavior. 

Data are constructed on an annual basis from 1974 to 1978 and are 
intended as extension of an existing data set developed at Rand (Sung, R-1665- 
ARPA) for the period 1957 to 1973. The coding rules and procedures of the 
preceding study are used in developing the present data set. The result 
is a set of data organized in a machine readable form and coded on computer 
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Biographical Information for each elite la coded on three carda. 

Card one containa Information pertaining to name, datea of birth and death, 
place of birth, generation, field-army affiliation by generation, commander 
or commissar by generation and by year, military region by generation and 
by year, functional affiliation by generation and by year, civilian or military 
background, and provlnce/ailltary region by year. Card two containa infor¬ 
mation pertaining to military, party, and governmental positions by year, 
and promotion and demotion (or both) frequencies by year. Card three, which 
pertains only to new elites included in our sample since 1974, contains 
information pertaining to year of birth, place of birth, generation, combat 
experience, award, civilian and military education. These background 
materials are available for the pre-1974 sample of elites in the Rand data 
set (Sung, R-1665-ARPA) and are not reproduced here. Table 6 displays 
the coding data layout for the computer cards. The coding rules used to 
register data maaerlcally on these cards are described below. 


3. Date of birth 

4. Place of birth (native place) 

I • Anhwei, 2 - Chekiang, 3 - Fukien, 4 • Heilungkiang, 5 " Honan, 

6 - Hopeh, 7 * Hunan, 8 ■ Hupeh, 9 ■ Inner Mongolia, 10 ■ Kansu, 

II - Kiangsl, 12 - Kiangsu, 13 ■ Kirin, 14 • Kwangsl, IS ■ Kwangtung, 

16 - Kweichow, 17 - Liaoning, 18 - Nlnghsla, 19 - Peking, 20 “ Shanghai, 
21 - Shansi, 22 - Shantung, 23 - Shensi, 24 • Slnklang, 25 - Sxechwan, 

26 - Tibet, 27 • Tientsin, 28 ■ Tslnghal, 29 • Yunnan, 30 • Taiwan 
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5. Generacion* 

GI - Pre-Hay 1928, G2 - June 1928-Novenber 1931. G3 - December 1931- 
December 1936. G4 - January 1937-Deceaber 1940. GS - January 1941- 
Auguac 1945, G6 - September 1945-Occober 1950. G7 - Noveober 1950- 
Sepcember 1954, G8 - October 1954-September 1959, G9 - October 1959 
-December 1963, CIO ■ January 1964-January 1967, Gil • February 1967- 
March 1969, C12 - April 1969-September 1973, G13 • October 1974- 

-September 1977, FE - Final Eatlaate 

6. Field-army (FA) affiliation by generation** 

1. - 1st FA, 2 - 2d FA, 3 - 3d FA, 4 - 4th FA. 5 - North China FA. 
China FA, 6 - the Center 

7. Commander (edr) or cooalssar (car) by generation*** 

1 • commissar, 2 - coomander, 3 • both 

8. Military or civilian 

1 - civilian, 2 ■ military 


*Generatlon affiliation la determined by the date of entry into 
the Army for military personnel, or into the CCP for civilians. In addition 
to the twelve generations identified in the Rand study we have added a 
thirteenth. The period of each of the generations represents more or less 
a cycle of crisis in the history of the CCP and the Red Army. Thus generation 
thirteenth Is Identified by the period bounded by the dates of the 10th 
and the 11th party congresses, roughly October, 1974 to September, 1977. 

••Field-army affiliation, comsnder or cotmsissar, military-region 
affiliation, and functional affiliation are data coded according to the 
generation period in which they occurred. 

•••Cosmander/coamlssar categories apply to military district, garrison, 
Rllltary regions, and service arms. In addition, the Director and Deputy 
Director of General Staff Department are classified as commander, whereas the 
Director and Deputy Director of General Political Department are classified 

as commissar. 


9. Military-region affiliation by generation 

1 ■ Slnklang, 2 • Kunming, 3 - Nanking, 4 - Canton, 
5 • Peking, 6 • Ch'engtu, 7 » Fuchou, 8 ■ Lanchou, 

9 - Shenyang, 10 - Tsinan, 11 ■ Wuhan, 12 - Center 
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10. Functional affiliation by generation* 

l - Armor, 2 - Artillery, 3 - Engineers, 4 - Infantry, 

S - Signal Corps, 6 • Railway Corps, 7 • General Political 
Department, 8 - 2d Artillery (Missiles), 9 - Public Security 
Force, 10 • Air Force, 11 - Navy, 12 - General Chief of Staff, 

13 - General Rear Service. 14 - Propaganda and Education (Including 
ministries: Culture, Education, Higher Education, Internal Affairs 
Public Health), 13 • Industry and CoamunlcatIon (Building Conatructlon, 
Building Materials, Chemical Induetry, 1st, 2d, 3d, 4th, 3th, 6th, 7th, 
and 8th Machinebuilding. Fuel 6 Chemical Industry, Commlcatlon, 
Geology, Light Industry, Testlle Metallurgical Industry, Petroleum 

Industry, Post and Telecoamunlcatlon, Railway, Water Conservation, 
Allocation of Materials), 16 - Finance and Trade (Commerce, Finance, 
Food, Foreign Trade, Economic with Foreign Countries, Planning 
Commission), 17 - Foreign Affairs (Foreign Affairs), 18 - Agriculture 
and Forestry (Agriculture, Forestry, Aquatic Products, Labor, State 
Farm and Land Reclamation), 19 • Political, Legal, Organization and 
Personnel 

11. Province or military-region affiliation by year** 

l - Anhwei, 2 ■ Chekiang, 3 • Fukien, 4 - Heilungkiang, 

3 ■ Honan, 6 ■ Hopeh, 7 - Hunan, 8 - Hupeh, 9 • Inner Mongolia, 

10 - Kansu, 11 • Klangsl, 12 • Klangsu, 13 - Kirin, 14 - Kwangsl, 

15 • Kwangtung, 16 - Kweichow, 17 • Liaoning, 18 - Nlnghsla, 

19 “ Peking, 20 “ Shanghai, 21 “ Shansi, 22 ” Shantung, 23 ■ Shensi, 

24 - Slnklang, 25 - Szechwan. 26 - Tibet, 27 - Tientsin, 28 - Tslnghal, 
29 “ Yunnan, 30 ■ Center MR, 31 “ Slnklang MR (Including 24), 32 " 
Kunming MR (Including 16. 29). 33 - Nanking MR (Including 1,2,11,20), 

34 - Canton MR (Including 16,15,7), 35 - Peking MR (Including 6,9, 
21,19), 36 - Ch’engtu MR (including 25,26), 37 - Fuchou MR (Including 
3,11), 38 - Lanchou MR (Including 10,18,23,28), 39 - Shenyang MR 
(including 4,13,17), 40 • Tsinan MR (Including 22,27), 41 - Wuhan MR 
(Including 5,8). 


*As distinct from the Rand study we collapse political-legal and 
organizational-personnel backgrounds into category 19, Instead of treating 
tham separately. We also Include In 19 vice premier, the minister and vice 
minister of National Defense Ministry, the Scientific and Technological 
commissar of the National Defense Ministry and the chief of the military 
academy. Moreover, categories 1 through 13 pertain national level People's 
Liberation Army organizations; categories 14 through 18 pertain to 
national level party and governmental organizations; category 19 also 
contains subnational level party and governmental conmisslons. 

•♦Provincial affiliations are coded at the lowest level for which 
data are available. The levels, considered In ascending order, are: 
province, military district, region, and center (the Peking district). 








12. Position level or rsnk, by year (1956-1973) 


Military (M) 


4 - National level 

41 “ Chairman, Military Affairs Committee (MAC) 

42 - Vice Chairman, Military Affairs Committee 

43 - Defense Minister 

44 • Vice Defense Minister 

45 “ Chief of service arm 

46 ■ Deputy Chief of service arm 

47 - Political Comissar (PC) of service arm 

48 - Deputy Political Commissar of service arm 

49 ■ Unidentified position 

3 - Regional level 

31 * Commander of military region 

32 • Political Commissar of military region 

33 • Deputy Commander of military region 

34 - Deputy Political Commissar of military region 

35 - Chief of Staff of military region 
30 ■ Unidentified position 

2 • District level 

21 • Commander of military district 

22 • Political Comissar of military district 

23 ■ Deputy Commander of military district 

24 - Deputy Political Commissar of military district 

Party (P) 

4 - National level 


41 • Politburo member 

42 “ Politburo alternate member 

43 ■ Pull member of Central Committee 

44 - Alternate member of Central Coasittee 


•Position Is identified at end of each year according to military 
region, functional afflliatlona, commander/commissar dlstrlnction, 
promotion/demotion/both categories. When no position is observed, no 
affiliation of any sort is assumed. In addition, we Include under 
unidentified position: (1) purged or demoted cases. (2) uncodable 
position, (3) missing data, (4) post-purge year. Uncodable position also 
includes cosmetic or nominal position to which no actual political authority 
is attached. 







106 


3 - Regional level 

31 • lac Secretary 

32 - Secretary 

2 • Provincial level 

21 • lac Secretary 

22 “ 2d Secretary 

23 ■ Secretary 

24 - Deputy Secretary 

25 " Standing aeaber 

Government (C) 

4 - National level 

41 - Chief and Deputy Chief of State 

42 - Premier 

43 - Vice Premier 

44 - Minister 

45 - Vice Minister 

46 - Director 

2 • Provincial level 

21 ■ Governor (Chairman, Revolutionary Committee) 

22 ■ Vice Governor (Vice Chairman. Revolutionary Committee) 

23 - Standing member, Revoluaclonary Committee 

98 - Dismissed and purged 

99 * Disappeared and attacked 
13. Promotion and Demotion* 


•Pr-’sotlon and demotion are coded according to movement between 
two successive years for the following cstegorles: position tj and 
position t., mllltsry region t. and military region t,. functional 
affiliations t. and functional affiliations t- province and province 
t 2 , conander/commissar and cotmunder/cosmlssar tj. 
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Card 3 

16. Year of birch 

17. Place of birth 

18. Year of entering party 

19. Coabat experience 

1 - Korea (1950-1953), 2 - Taiwan Straite (1955, 1958, 1962), 

3 - Sino-Indian border (1962), 4 - North Vietnam (1965-1971), 

5 • Slno-Soviet border (1965- ) 

20. Awards 

1 - August 1 medal, 2 - Independent Liberty medal, 3 ■ Liberation aedal 

21. Civilian education 

1 - grade school, 2 - high school, 3 • college, 4 - study 

in USSR, 5 - study in Germany, 6 - study in Japan, 7 - study in 

France 

22. Military education 

1 “ basic school, 2 - anti-Japanese college, 3 ■ Nanking 
Staff College, 4 - Peking War College, 5 - Soviet advisers 
(including study in Soviet ailltary schools) 

According to these categories we group the elites into factions. 

The total nuaber of elites in each category represents the base score froa 
which the crisis indicator is constructed. Assuming that a large nuaber of 
deaotlons signifies abnormal strain on the political system and thus the 
occurence of crisis, we count the total frequency of deaotlons in each 
elite faction. Coaparlng the demotion frequencies for competing factions 
gives us a further measure of inter-factional conflict. 

Deaotlon is taken here to aean not only aoveaent froa a higher 
position to a lower one but also purge, dismissal, and disappearance. To 
coelute deaotlon we coapare the lowest position held by an individual for two 
consecutive years. If the position held in the earlier yaar is higher than 
that of the latter than e deaotlon is assumed to have occurred. (In contrast. 
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a movement In Ch« reverse direction signifies e promotion). Demotion 
frequencies ere counted for elite factions, within the party, the governmental 
and the military hierarchies. A specific number of such positions have been 
charted in the previous Rand study (Sung, R-1665-ARPA) through which 
demotion or promotion patterns may be traced. As the procedure is clearly 
described in the preceding study no further attempt will be made here to 
explicate the counting rules. All further explanations are given in the 
previous study. 

The policy indicators are developed from aggregate economic data. 

With the exception of the agricultural and industrial output measures most 
data are estimated in constant dollars. These Include Cross National Product, 
Import, export and defense expenditures. Since annual agricultural and 
Industrial output measures are difficult to obtain and not available for many 
years, only estimates of performance levels provided by the Central 
Intelligence Agency are used. These measures are derived using the 1957 
production figure as the base year to which a value of 100 Is assigned. 

Against this base value the output performance levels of all other years are 
measured. Thus, the performance indicator measures for other years may be 
higher or lower than 100 depending on their levels of output relative to 
that of the base year. 

The sources for the above data are listed in the following references. 







Bibliography of Oaca Sources 


Appearance* and Activitiea of Leaning Personalities of the People's 
Republic of China . 1 January - 31 December 1976, (May, 1977) 

Central Intelligence Agency, Washington, D.C. 

Appearances and Activities of Personalities of the People's 

Republic of China , 1 January - 31 December 1974, (May, 1977) or 
FF12059, Central Intelligence Agency, A(CR) 73-10, Washington, D.C. 

Appearances and Activities of Lading Personalltltes of the People's 

Republic of China, 1 January - 31 December 1973. (March, 1976) Central 
Intelligence Agency IR76-10980, Washington, D.C. 

Biographic Service . Union Research Institute, Hong Kong. 

"China: Cross Value of Industrial Output, 1963-1977." Central Intelligence 
Agency (1978). Doc«ent Expediting (DOCEX) Project, Library of 
Congress, Washington. D.C. 20340. 

"China: .International Trade 1977-78." Central Intelligence Agency (1978). 

Document Expediting (DOCEX) Project, Library of Congress, Washington, D.C. 
20540. 

"China: Economic Indicators." Central Intelligence Agency (1978). Document 

Expediting (DOCEX) Project, Library of Congress, Washington, D.C. 

20540. 

China Directory . Raido Press Inc., Tokyo, 1977. 

China Directory . Radio Press Inc., Tokyo, 1979. 

Chinese Communist Affairs Monthly . Center for International Studies, National 
Chengchl: University, Taipei (January 1973 - June 1979). 

Chinese Comaunlst Who’s Who , Vola. 1 and 2, Institute of International 
Relations. Republic of China, Taipei, 1970. 

Chinese Co nan uniat Who’s Who . Revised edition. Center for International Studies 
National Chengchl University, Taipei, 1978. 

Chung Rung Nlen Pao (1973). Institute for the Study of Chines* Comammlst 
Problems, Taipei, 1973. 

Chung Rung Nlen Pao (1974). Institute for the Study of Chinese Conunlst 
Problems, Taipei, 1974. 

Chung Rung Nlen Pao (1975). Institute for the Study of Chines* Communist 
Problems, Taipei, 1975 

Chung Rung Nlen Pao (1976). Institute for the Study of Chines* Communist 
Problems, Taipei, 1976. 
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Chuns Kua| Nlen Pao (1977). Institute for Che Study of Chinese Communist 
Problems, Taipei, 1977. 

Chung Kung Mien Pso (1978). Institute for the Study of Chinese Communist 
Problems, Taipei, 1978. 

Directory of Officials of the People's Republic of China . Central 
Intelligence Agency, A73-35, Washington, D.C. January 1974. 

Directory of Officials of the People's Republic of China . Central 
Intelligence Agency. A(CR)75-16, Washington, D.C., April 1975. 

Directory of Officials of the People's Republic of China . Central 

Intelligence Agency. CR77-15208, Washington, D.C. October 1977. 

Directory of Officials of the People's Republic of China . Central 
Intelligence Agency, CR78-11373, Washington, D.C. April 1978. 

Directory of Officials of the People's Republic of China . Central 

Intelligence Agency, CR78-16506, Washington, D.C. November, 1978. 

Hierarchies of the People's Republic of China . Union Research Institute 
Hong Kong March 1975. 

Huang Chen-Lsia, Mao's Generals . Research Institute of Contemporary History. 
Hong Kong, 1968. 

Issues and Studies . Center for International Studies. National Chengchl 
University, Taipei (Jan. 1973 - June 1979) 

Klein, Donald and Anne Clark. Biographic Dictionary of Chinese Communism : 
1921-1965 Vols. 1 and 2, Harvard University Press, Cambridge, Mass., 
1971. 

Malcolm Lamb, Directory of Chinese Officials and Organizations, 1968-1978 . 
Contemporary China Center. Research School of Pacific Studies 
Australian National University, Canberra, Australia, 1978. 

Ting Wang , Chung-Kung Wen-ke Yun-tung Chung te Tsu-Chlh YU Jen-shlh Wen-tl . 
1965-1970 . Contemporary China Research Institute, Hong Kong, 1970 


Who's Who In Communist China . Union Research Institute, Hong Kong, 1966. 

"World Military Expenditures and Arms Transfers, 1967-1976," United States 
Arms Control and Disarmament Agency (1978) U.S. Government Printing 
Office 1978. 
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Statistical Model for Hypothesis Testing 

Data analysis tests the relationship between five crisis and four policy 
indicators. These indicator are defined in sections 4.4.1 and 4.4.2. 
Consistent with hypothesis H., we test for causal effect of crisis on policies, 
as well as the reverse feedback effects. The policy indicators are examined 
one at a time in relation to all the crisis indicators. The reverse feedback 
effects are, in turn, examined. Using the autoregressive-moving average 
vector model our estimates of relationships for the defense expenditure 
indicator are: 

T-k 

(defense expenditure/GNP). ■ Z (Commander/Commissar) 

1 t-T llC C 


I (field army i/field armv J) 

t-T 1 


♦ Z (function 1/funtlon J) 

t-T 1JC C 


(generation i/generation J) 


(civilian/military) 
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(commander/commissar) T 
(field army 1/fleld army J)^ 
(function 1/funtlonj)^ 
(generation 1/generation J)^ 


(civillan/mllitaryj^ 


T-k 

t 

t-T 

T-k 

t 

t-T 


*21t 

*31t 


T-k 

I 

t-T 


*4lt 


T-k 

r 

t-T 


*5lt 


T-k 

I 

t-T 


*61t 


(defense expenditure/GNP) c + N2 C 
(defense expenditure/GNP> t + N^ t 
(defense expendlture/GNP> t + N^ t 
(defense expendlcure/GNP)^ ♦ 
(defense expendltures/GNP) c + 


where Is a moving average process of the form 
T-k 

N lt “ t f T 9 ljt “it ,1-1.6; j-l.5 


The first equation above defines the effects of the crisis Indicators on a 
policy Indicator, whereas the remaining five equations define the feedback 
effects of the policy on crises. These relationships are repeated for three 
other policy Indicators: lndustrlal/agrlcultural outputs, export/GNP, lmport/GNP. 

Alternatively, we have also entertained models In which possible 
Interdependence of the policy indicators are considered. Thus, a policy Is 
assumed not only to be affected by crises but also by other policies. The 
model Is a system of autoregresslve-movlng average vector functions. In this 
analysis only one crisis Indicator is considered at a time In conjunction with 
all of the policy Indicators as the following vector expression illustrates. 







(defense expendlture/GOT), 
(Industry/agriculture) 
(laport/C»P) T T 

(export/GNP) T 
(commander/commlssar)_ 


Some results of estimations for the 
51. 



(defense expenditure/GNP)^_ j 
(Industry /agriculture).^ 

(import/GOT)^ 

(export/GNP) t _j [ 

(commander/commissarj 



(defense expenditure/GNP)^ 
(lnduetry/agrlculture) T _ k 
(Import/GOT)_ . 


(commander/commissar)_ 



models are displayed on tables 5a to 
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